{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "# from __future__ import print_function\n",
    "import pandas as pd\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import statsmodels.api as sm\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0    597816.0\n",
       "1    583104.0\n",
       "2    572465.0\n",
       "3    561279.0\n",
       "4    551589.0\n",
       "Name: number, dtype: float64"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "discfile = '..\\\\dataset\\\\data_test.csv'\n",
    "data = pd.read_csv(discfile)\n",
    "data=data['number']\n",
    "data.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.axes._subplots.AxesSubplot at 0x502bb38>"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZEAAAD8CAYAAAC2PJlnAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJztvXl4Y1d9uP8eyZIlWZb3ZcaefV9CJjOTPYQsZIVmUiAtFMhQ0oZCCqXQFmj7NA9LW/g+FNrwawOBJCQ0ZWkgTUgDYTIhoZB1JpnJMpOZ8axeZrxbtixZ6/n9ce+VZY9sS7JkXeme93n8SDo6995zZel+7mcXUkoUCoVCocgFW7EXoFAoFIrSRQkRhUKhUOSMEiIKhUKhyBklRBQKhUKRM0qIKBQKhSJnlBBRKBQKRc4oIaJQKBSKnFFCRKFQKBQ5o4SIQqFQKHKmotgLyDeNjY1y+fLlxV6GQqFQlBR79+4dkFI2Zbtd2QmR5cuXs2fPnmIvQ6FQKEoKIcTJXLZT5iyFQqFQ5IwSIgqFQqHIGSVEFAqFQpEzSogoFAqFImeUEFEoFApFzigholAoFIqcUUJEoVAoFDlTdnkiCoVCUU785nA/zb5KBgMRLlnVgBCi2EuaQkaaiBDiL4QQbwgh3hRCfFofqxdC7BJCHNEf6/RxIYS4SwjRIYR4TQixNWU/O/X5R4QQO1PGtwkhXte3uUvon9JMx1AoFAorEIsnuPW+l7j+X/+PD37vRZ451F/sJZ3FnEJECLEZ+FPgAuBc4N1CiDXA54HdUso1wG79NcANwBr973bgbn0/9cCdwIX6vu5MEQp363ON7a7Xx2c6hkKhUJQ9B06PTnl9cnC8SCuZmUw0kQ3AC1LKoJQyBjwL/D6wA3hAn/MAcLP+fAfwoNR4AagVQiwCrgN2SSmHpJTDwC7gev09n5TyeSmlBB6ctq90x1AoFIqy5/Vu/5TXR/tLU4i8AVwuhGgQQniAG4ElQIuU8jSA/tisz28DOlO279LHZhvvSjPOLMdQKBSKsmd4PALAvn+4hs1tPo72B4q8orOZ07EupTwohPgamuYQAPYDsVk2Sef1kTmMZ4wQ4nY0cxhLly7NZlOFQqEwLf5QFLfDTq3HybL6Kg6eGZ17owUmI8e6lPJeKeVWKeXlwBBwBOjVTVHoj3369C40TcWgHeiZY7w9zTizHGP6+u6RUm6XUm5vasq6krFCoVCYkpFglFqPA4Cm6kr6R8NFXtHZZBqd1aw/LgXeA/wQeAwwIqx2Ao/qzx8DbtWjtC4C/Lop6kngWiFEne5QvxZ4Un9vTAhxkR6Vdeu0faU7hkKhUJQ9I6EoNW5NiLT4XIyFYwQjsxmCFp5M80R+KoRoAKLAHVLKYSHEV4GfCCFuA04Bt+hzn0Dzm3QAQeCPAaSUQ0KILwMv6/O+JKUc0p9/HPg+4AZ+of8BzHQMhUKhKHv8KUKkuboSgL7RMMsbzZPil9FKpJRvTzM2CFydZlwCd8ywn/uA+9KM7wE2Z3oMhUJRukgpTZcwZ1b8wSjLGjwANPt0ITIWZnljVTGXNQVV9kShUCwYP3ulixVfeILdB3uLvZSSYCQUSfpEWnwuAPrGJoq5pLNQQkShUCwYP3zpFAD7u/xzzFS82eOndzTM2pZqYNKc1Wsy57oSIgqFYkGQUnLozBgAZ/yhIq/G/BglTm7ZpgW11rgdOCtsShNRKBTWpGs4xOiEFll02m+uC6EZGRqPUOW0U6Obs4QQNJswzFcJEYVCsSC82aMlyi2qcXFGCZE50XJEnFPGmqsr6VWaiEKhsCIHevzYbYJrN7ZwcjBI13Cw2EsyNSPBSae6QVN1JQNjkSKtKD1KiCgU8+CLP3+Tj37/ZcbD5koAMyNvnRljRWMVt79jFdFEgp/u7S72kkzNcBohUuN24A9Fi7Si9JgnY0WhKDGklNz/uxOAdoHctky1u5mNzuEQy+o9tNW6WVzjNmUxQTMxEoqyqNY9ZcyMQkRpIgpFjpwZnbRNDwbM5ew0I93DQdrqtIviyqYqjg+Yr6y5mRgJRqlLo4mEonEisUSRVnU2SogoFDlyLKW3w9C4uezUZmNsIsroRIw2/c56VZOXY/0BtAIXiukkElLzibinOtaNEihm0kaUEFFMoWs4SCxunrscM5N6Jz2ohMis9IxoWtvi2klNZDwSp29MaXDpGAvHSEjO8on4lBBRmJmu4SCXfe3X7Lz/pWIvpSToGQlRYRO4HDaO9I4VezmmZnBcExaNXi3reoVe+0n5RdIzEtRuSqaH+CpNRGFqfvnGGQB+1zFY5JWUBt0jIRbVupiIJviffT281jVS7CWZltGQFr3mc2uxPCubvMBUk6BikpGgJiTS+UQARpUQUZiRF49rlfltqsBqRvSMhFhc46ZJr2n06iklRGZidEK76Plc2kVwkc+Fs8LGqSGVK5KOYaWJKEqRA3pGcUJCPKEcnnPRPRyirc7Nb/76Spx2Gx19yjQzE8ads2HTt9kETd5K+pVPJC2GJpIuTwSUEFGYEH8oSvdIiOpKzdwQVc71WZmIxjk9OsHyhircTjub23wcVn6RGRmdiCEEye8X6O1elRBJi+ETqZumiSjHusK0dOpmhbWtWtlpVSBvdo4PjCPlpIO4rc5D76j6zGZiNBTFW1mBLcVW2qyEyIwMBw3z39R8cIfdRpXTroSIwnwYQmOJngx25defYSIaL+aSTI3hEF7ZpAmR5upKekfDKu9hBsYmYkl/iEFTdSX9KkkzLf5QlGpXBRX2sy/RZstaV0JEAUxmXy+p90yOKW1kRnpGtH4YxufVXF1JKBonoGpopWV0Ipo0xRg0VVcyNB5RptM0pPZWn47P7Uj6TMyAEiIKQGsSZLeJZAtOmLxQKs5mbCKKEOB1auaGydal6s46HUPjEWqnXRTrqzR7v5kuiGbBH4qe5VQ3qHE7VIivwnyc8Ydprq6ksmLyK9GlhMiMjIVjeJ2TNv7J1qVKe0vHycEgS1O0XEjJeZgwzwXRLIwEIzNqIrUeZc5SmJCh8TANXieOFBts97ASIjMRmIjhTXF6Nvs0IaIcxWczHo4xEAiztGGqEDF8JGa6qzYLs5mz6qsqkxUAzIASIgpg8ktbYZ+MnlElKWYmEI7hrUwVIro5y2StS82AkVC4bLoQ0bPXzXRXbRb8oRg104ovGiyqcTEQiBCOmSPwRQkRBTApRFI1kYOnR4u4InMTCE/VRKorK3A5bMqclYbTfk2jXZymNwaQ7Luu0JBS4g/NbM5qrTHXDYsSIgrAuPNx4EjRRI4PjKsw3xkYm5iqiQghaK52Kcd6GgxNY7pj3TBnKU1kKqFonGhczihEFulCxCy5XEqIKJBSMhrSQjArbJNfiYRUNv6ZCIRjVE9LBGvxVdI3Zo4ftpnw69FX0y+KPhMWEzQDSaE7Q3TWpBAxh89SCREFE9EEkXjiLJ8IqPDLmQhM00QATRMxiYnBTPiTFXynXhRdDjvOCpuKzprGyAxC16ChSgviGAyYo4eNEiKK5J3PdJ8IwFDQHF9UszE2EaVqmhBpqq5U5qw0+PWSJ9O/W6D5ksaUT2QKM5n/DAxhPGISDU4JEUWy7HQ6ITKihMhZ+INRxiPxpFnBoMXnIhCOMa6y1qcwW7iq22lXfrdpGJrIdM3NwG4T+FwVpjEDKiGi4I1uPwBrW6qpmNZMZFi1fT2LowNa6PPKRu+UcSPhUGkjU/GHzi55YuByKCEyndHQ7OYsgBqPwzQ3eEqIKHjp+BB1Hgerm7xTNBEhJquJKiY5qvcNMYovGhgJh30qzHcKo6EoNe6KtO+5HXZCESVEUpnLsQ5Q63aWljlLCPGXQog3hRBvCCF+KIRwCSFWCCFeFEIcEUL8WAjh1OdW6q879PeXp+znC/r4ISHEdSnj1+tjHUKIz6eMpz2GIr+cGgqyprkam01Mcaz7XOa52zETncMhhJharBI0xzooTWQ6s5qzHHZCShOZwkgogt0mzgrcSMVMlXznFCJCiDbgU8B2KeVmwA68H/ga8E0p5RpgGLhN3+Q2YFhKuRr4pj4PIcRGfbtNwPXAfwgh7EIIO/DvwA3ARuAD+lxmOYYij4xOxJLmBkdKiG+dx8GQ0kTOYnQGR7FxoVSO4qnMJkRcTjuhqKrim4o/FMXnqkCImftU13gcydDpYpOpOasCcAshKgAPcBq4CnhYf/8B4Gb9+Q79Nfr7Vwvt09gB/EhKGZZSHgc6gAv0vw4p5TEpZQT4EbBD32amYyjyiJYjot31pGoitR6n0kTSMDoRPas3BpDMYA+EzfHjNguzayI2wkoTmYI/FDurt/p0aktJE5FSdgNfB06hCQ8/sBcYkVIat1xdQJv+vA3o1LeN6fMbUsenbTPTeMMsx1DkkdSLYurddZ3HkYzcUkwyGoqldRR7HHaE0HJIFBrhWJxQNK7MWVkwEozMGIhgUON2MBKKmqIJWibmrDo0LWIFsBioQjM9Tcc4m3Q6mMzjeLo13i6E2COE2NPf359uimIGEglJIBxLtuE0yp7YBNRVORkeN8fdjpkYnYiela0OYLMJqpwVBMLqomjgnyPSyKUc62cxOovmZlDrcRDXf7vFJhNz1juB41LKfillFPgZcAlQq5u3ANqBHv15F7AEQH+/BhhKHZ+2zUzjA7McYwpSynuklNullNubmpoyOCWFQSASQ8rJmHSjHacQgjplzkpLulavBt7KCmXOSsEIV50txFdpIlPxh6IzJhoaGELGDCatTITIKeAiIYRH91NcDRwAfg28T5+zE3hUf/6Y/hr9/aelpnM9Brxfj95aAawBXgJeBtbokVhONOf7Y/o2Mx1DkSeSP3L9omjkidiEZs4aj8RNU3LaLKT6kKZTVWlnXGkiSebSRFSy4dmMZKCJGGXizVCWKBOfyItozu1XgNf1be4BPgd8RgjRgea/uFff5F6gQR//DPB5fT9vAj9BE0C/BO6QUsZ1n8efA08CB4Gf6HOZ5RiKPDGq1zUyzDNOXRP5/A0bks49M3xRzcRMjnUAr8vBmAlMDGZhTiHisBONS2KqzzqgmZczNWeBOYpXzhyInIKU8k7gzmnDx9Aiq6bPnQBumWE//wj8Y5rxJ4An0oynPYYifxjF7wxzg80mOPHVdwHw2H7Nejg2EaPFV5z1mY3pPqTpeCvtquxJCpkIEYCJWAJvmtpaViMQiZGQsycawuTnaYaEQ/VfszgjKXWzpuPRf+DBiLooGkz3IU3HW1mhorNSmKkMvIHLqb5jqfjnqJtlYAgZM1gJlBCxOMaXsK7q7Lh0T/IHrmzWBtN9SNPxVjoYU6XNk8xUBt6gRa831jlkjt4YxWYuzc2g1BzrijLGqI2VLhrErQsRFYI5yXQf0nQavU4GAhFTxO+bAX8oSpXTnrYMPMCWpbUAvHpqeCGXZVpGZvk9puJ22HHabYyEih89qYSIxRkJRXDabUmtIxWPU7tQKk1kkuk+pOm0+FxE4glVuFJntmx10OqNLa5xJStJW52kJjKHT0QIQY3HYQrHuhIiFmdkPEqNx5G2To9H2avPwqiLNZM5q1XvMXLGJP2vi81sZeANFtW6VdFKnUzNWcYc5RNRFJ2RUIS6Ge563MonchaTyXPpzVktPk2I9Kpe60Bm2ddN3kr6lRABSJqnat1zFyyvVUJEYQaGg9EZv7DKsX42hjmregZNpEXvKdKrNBFgbnMWaG2F+wNKiID2eTntNlyOuS/NZikHr4SIxfEHozPGpLsqDMe6MmcZzOVYN3qKnFGNqYDMhchIMKoqI2BUQ0hvXp5OjUcJEYUJGA5GZhQiNpvA7bArTSSFsYkonlmijZwVNhq9TnqVEAEyEyKNXk17GwwUP9Ko2IzMclM3nVq3OWrbKSFiYaSUjISi1M3Su6Cq0k5Q1TZKMlMF31Saq130jirzTCSWmLUMvEGjV/v+DSiTVkZC18DrqmA8EieRKG44uRIiFiYUjROJJWZtgON22gmqMh5JRkMzV/A1aK1xqegsMg9XrdcTXVVYdHZCxCgZE44Vt+6YEiIWJpnYNMuPvM7jZHC8+CqzWRgLzx2y2uJzKXMWmYerGtUSzGCaKTYjwbnLwBu4ded7sasgKyFiYYyuhTOF+AIsrfdwaii4UEsyPZomMpc5q5LB8QhRi1em9YcyqwNlmFOH1M1K0rGeCS5dEyl2PxYlRCzMZHG8mc1Zyxo8dA+HVKluHc0nMvuPvFkP87W6jd8fmrm4Zyo1bgdCwLDFhUgsnmAsHMvcnOVUQkRRZDIxNyyrryKWkPSMKPMMzN6QysAI8+2zuHM9U3OW3SaodTsYsrg5a1SvhpBpdFZSEyly9KQSIhbGaJ40W7TR4lo3AKf9qsqqlHLW1rgGzXplWquX8jB6zVdXzt22qM7jtLxjPZuSJ5DSi0VpIopiYTRP8s7yI2+qNkwz1r5LBM1sEEvIOW3Whjmrz+KlT4wkVXea4p7TqatyWt6cZfweqzIQujD5uU5EVXSWokgYzZNm+9IaMfz9Fr8gwmTxxdmELkwmz1m9HlQool3cjGrQs1HncVresW5oFIaGMRdu5VhXFJtAJIazwoazYuavQZ3Hid0mlCbCpBCZK9nQYbdRX+W0vDkrGNW+X3bb3CU86qscyWhBq2IIg3RtGdKhorMURScwEZvTXm2zCeqrnJaPNAIIZOBDMmiurrS8Yz0UiWd8V234RKzczMtwkLsy/MyMIo0TyrGuKBbj4VhG9ldVqlvDaHvrrZzb8dlUXWl5E2AoEs/4rrquykkklrB0nTZDo8jEhwTKnKUwAYFwbE77PkBjdaXSRJj0IWWmibiUOSsaz/iCWK8SDrP3iag8EUWxGZuI4c3ggmj0Dbc6YxlEsxk0+zTtrdjF8YrJRDbmrGT9LOt+zwxzVqafmdGqodgVJZQQsTDjkcw0kSav1jTIyvZqyE4Tqfc4iSUkAQv3YglmY87SE+ysnCsS0kN1M9XebDbBTecu5scvdybDg4uBEiIWJjCRoRCpriQSSyQzaq1KIIs4fiPreGTcuhdFzZyVWc5DUhOxsDnLMEtVzhItOZ3L1jQST8iimgGVELEwgXA8owuikfdgdb9IIBzD5bDN2JAqFaO8vtEz24po5qzMLjHKJ6L5RNwOe0ZdDQ3M8LkpIWJhAuG5GyxBihCxuKN4bCKWUWQWTJpnRixsnglGYxklGoJW6dcmlE8kU1OWQb1XCRFFkYjFE0xEExlGZ+lZ60oTyUjowqQ5S10UM7so2m2CGre1Ew5D0cwDEQyUJqIoGuN6cbxM80RAaSKBiWhGQhcmy+sbRfWshmGnb6iauc3AdLT6Wdb8vEATuq4MzX8GhiZSTOGrhIhFMaKGMq2wqkqf5KaJWNWcNRAIk5Bal8dMsXL9rIFAmKff6ku2EciU6soKHHZR1O6jSohYFCNcNZM8EaP0idWz1scyjGYDrX6Wt7LCsuYZo8d8NkLE47QXPXGuWBw8PUooGuejl63IajshhFYyxsxCRAixTgixL+VvVAjxaSFEvRBilxDiiP5Yp88XQoi7hBAdQojXhBBbU/a1U59/RAixM2V8mxDidX2bu4QenjDTMRTzJ5twVdBMWlaPzso0OdOg1uNIdo+0GkaP+dYshEhlhY1IzJodNA2NdVmDJ+tt66uc5tZEpJSHpJRbpJRbgG1AEHgE+DywW0q5Btitvwa4AVij/90O3A2aQADuBC4ELgDuTBEKd+tzje2u18dnOoZingSyyL4GaPA6GbCoqcEgEJ67YGUqtR4HIxb1ifTqWmuL3lslE5wVNiIWbcNs+M5qM2xIlYrpNZFpXA0clVKeBHYAD+jjDwA36893AA9KjReAWiHEIuA6YJeUckhKOQzsAq7X3/NJKZ+XWkr0g9P2le4YinkSyLA3hkGN28GYRS+IoHU1DISz1ETcTsuaswZ1rbU+C8d6ZYWdcMya5ixDiMzV8Cwd9V5nUVsLZytE3g/8UH/eIqU8DaA/NuvjbUBnyjZd+ths411pxmc7hmKeGBe3TPs5+9wORiesK0QmogniCZlxnghY25w1PB6hxu2gIoPETAOn3brmLH8oisthy7gMfCr1RQ5IyPg/LIRwAjcB/z3X1DRjMofxjBFC3C6E2COE2NPf35/NppbljH8Cu00kEwnnotpVYemyJ2NhvQx8lj4Rq2oiw8FoMuEyU5yW9olEMu6tPp36Kif+UJRYkUyB2WgiNwCvSCl79de9uikK/bFPH+8ClqRs1w70zDHenmZ8tmNMQUp5j5Ryu5Rye1NTUxanZF16/CFaqisz6joH4HM5iMQSyXLVVmM0pAlQXxZCpM6j/bitWMl3OBhJ1sPKlMoKG2GLChF/KEqtO7vPy6C+yomUFM3/lo0Q+QCTpiyAxwAjwmon8GjK+K16lNZFgF83RT0JXCuEqNMd6tcCT+rvjQkhLtKjsm6dtq90x1DMkzP+CVprMo+cMWy1VjVpGTb+TDU30PxICTlZQt5KDAcj1HmyuyhaWxOJ5qyJGLlLgSJZCjISIkIID3AN8LOU4a8C1wghjujvfVUffwI4BnQA3wU+ASClHAK+DLys/31JHwP4OPA9fZujwC/mOIZinpzxT7Coxp3xfOMOfMyiJi0j0bLBm/mFMVmE0YImreHxaNZCpLLCTiwhiVtQcxsaj2T13UrFCNMPFOlmJSPdXEoZBBqmjQ2iRWtNnyuBO2bYz33AfWnG9wCb04ynPYZi/vSPhbl8beZ31T6XrolYNEJrcFzTRBqqMv/MUoswLmuYY3KZoWki2ftEACKxRNaFCEudwfEIF2Rp/jMwIiyL1VNEZaxbkHAszlg4llVdI59b+6Ja1bk+MBZGiOxCVpOlTywmeCeicYKReNY+kVQhYiVi8QTDwQgNWZhKUzE0kfEiNUBTQsSCGEXu6rNQn62uiQyMR6jXa4hlilXNWUZEWvbmLO1yFI5bK3hjOBhFSq0NdS54KzWtLRAuzuemhIgFmTTNZKOJaELEqj6RoUAkKy0EJrOPrVaEMXmTUpWbOeuif9qd9zWZmVxMpalUKXOWYqExfuTZ3CkaESBWjc4aCUUyTsw0qLGqEEkmsuamiVjNrz6kB21ke5NioISIYsFJ3vlkoT67HXYqbMKy5ix/KJZ1CGaF3Ua1y3qVfI3zzfaimNpb3Eq5NUZNulzNWVXO4kZnKSFiQYwSCfVZqM9CCHxuh2XNWaOhaE51jWo9Dss1pjKKAWaruTlThIiVghGMHKRcHet2m8DtsCtNRLFw9I+FcdhF1hVDfa4Ky5qz/KHcksHqPNYrwjgcNCrSZndnnRqVZaW2A4OBCDaRWwVfg6rKCuVYVywcfWNhGr2V2LKINAKodjksac6KxRMEwtmbs0Dzi1jNJzI0HqG6smKKZpEJqRqblRqgDY5rQRvZ/h5TWVrv5s0efx5XlTlKiFiQ/rEwTdXZq84+d4UlzVnGOeciRGr1+llWYiQYoTbLyCyAK9ZNFum2liYSzjkyy+DaTa281uWneySUp1VljhIiFqR/LExTDvZXn8uaVWmTvR5cuZizrPeZDQWj1GcZmQVaK939d14LWE8TybXkicGNmxfxF1evwWHPXZvJFSVELEh/IExzFh3nDJY1VNE5FCJqse5zhpM3W0cxaHZuq1XyHQlGsg7vNfC5NDNYv9U0kRyd6gZLGzz85TVraa7OvKhqvlBCxGLEE5LBQG6ayNoWL5F4gpOD4wVYmXmZT/RMrUcr020lM+DQePaJmQZCCJq8ldbTRHL8vMyAEiIWY3A8TEKSk09kbUs1AIfOBPK9LFMzkCwDn/0PfWm9B4AXjg/mdU1mZiQYzUlrM2isto4QCcfijE1kV8fObCghYjGMH2cuQmR1sxebgEO9Y/lelqkxysBn00vE4Ip1TTRVV/K/r53O97JMSTgWJxCOZV03K5Umb2XyMy93jJyt+ZqziokSIhZjPkLE5bCzvKGKI5YTImG8lRU59b+usNtor3Nbxrk+qF/8c/l+GTRVOy2jiQzm0KfGbCghYjH6DCHizc0Bt6bFa0lNJNeSFKA5162SK5K8SZnHnXWTt5Kh8bAlmlMNzrPkiRlQQsRizEcTAWir9dDrn8jnkkzP8DwcxWCtXJH5fr+MbRNy0tRTzpzxa3kdud7UmQElRCxG51CQRq8z585xdR4H45E44Zh1ej6MR2LJSqm5oGWtl/8FEUiG5s5HiBi+JyuYtI4PBHHYBYtrlRBRlAgnB4Ms0SOGcqG2ymi0ZI07a4BQJI47B3+IQY3bwehEzBLmGePCPx8bvyGArJC1fnwgwNJ6DxX20r0Ul+7KFTlxaijIsnkIESMT2SqOYoBQND6vnt9GuKsV6o71j4Wp9TiorMj98zKEiBU0kRMDQVY0eou9jHmhhIiFmIjGOe0PsbShKud91OkXRKOxlRUIRuJ45iFEjJpbVhC8uZbUScUwZ5W7JpJISI4PjrOyKfffoxlQQsRCHDozRkLCxkXVOe/Din3DNXNW7j4Rwyk/aAFHcX8gt+KeqVRVVuBx2steE+nxh4jEEqxoVEJEUSK82TMKwKbFNTnvw7ggDllEiEgpdXNW7j+V1hrNaXrGAlFtuVaInk6jt7Ls62cdH9DKBykhoigZ9pwYos7joL3OnfM+DPu+VRzrkXiCeELiceauiSzyaZ9376g1hEgumf3TaaquLHtzlhIiipIinpD8+lAfV6xrRojcy0W7HHbcDnuyBWq5E4pooczzic7yuStwOWxlr4mMh2OEovG8aCJWKMJ4fGAcj9NOcx4+r2KihIhFeOXUMMPBKFdvaJ578hzUVzmTLVDLnaAuRObjWBdC0OpzcabMNZHTeuJcq2/+OQ+N1c6yr591fGCcFY1V87qpMwNKiFiEpw72UmETXL62ad77qrVQo6VQVNdE5iFEQDPPDJb5RbFzSBMiS+pzN5caNHldDI1Hyrp3jSFESh0lRCzC7oN9XLiyPqfufNOp8zgtI0R+8bpWfXc+5izQotrK/TPrHA4CsKQu9zwkg8ZqPaKtTAVvOBancyjISiVEFKVAR98YHX0Brl7fkpf91XqsU1Dw6786DICjYn4/lToLfGZdwyGcFbb8ONbLPFfkxECQhIRVzaWdaAhKiFiCB5+Ze+1qAAAgAElEQVQ/icMuuGnL4rzsr77KaYnieFJOlik5py33sGiwiCYyFKS9zo3NNn8bf7lnrR/t1xq7rWpSQkRhco71B/ivF0/x3q3teblDBGivc+MPRcs+QmtUb2n7dzdumPdnV+txEI4lktFe5UjncDAvpixIKcJYpprI0T4lRBQlQCye4O8eeQOXw85nr12Xt/2ua/UB8NaZ8u4rYkQbLa6dv6O4zgI1x7qGQ/PKQUrF0ET6yjSiraM/QFute94BG2YgIyEihKgVQjwshHhLCHFQCHGxEKJeCLFLCHFEf6zT5wohxF1CiA4hxGtCiK0p+9mpzz8ihNiZMr5NCPG6vs1dQo95m+kYisz4wQsnef7YIP/wexvzErtvsKFVK5vy1pnRvO3TjPSMaEJkUR7KdCdrjpWpEBmbiDISjM6rQnQqLoedRq+Tbv1/UG4c7Q+UhT8EMtdE/g34pZRyPXAucBD4PLBbSrkG2K2/BrgBWKP/3Q7cDZpAAO4ELgQuAO5MEQp363ON7a7Xx2c6hmIO/MEo/7b7CJeubuCWbe153XdTdSVOu63s8x56RrTzW1yTP02kXH1JRnhvvjQRgLZaN13D5SdEEgnJ0b5xVpeBKQsyECJCCB9wOXAvgJQyIqUcAXYAD+jTHgBu1p/vAB6UGi8AtUKIRcB1wC4p5ZCUchjYBVyvv+eTUj4vNU/mg9P2le4YilmQUvI3P93P2ESMv71xQ96TmYQQWoRWmVfy7RkJUWETedHimvUEvL7R8rTxd+UxvNegrc5dlprI6dEJQtE4q5pLP7wXMtNEVgL9wP1CiFeFEN8TQlQBLVLK0wD6o5EK3QZ0pmzfpY/NNt6VZpxZjjEFIcTtQog9Qog9/f39GZzS2fzbU0f4yuMHctrWbNy1u4Mn3+zlCzesn1exxdnQstbL867a4LR/ghafC3seoo2M0ha9Y+WpvXUOG4mGeRQitW56RkJTouTKgXJyqkNmQqQC2ArcLaU8DxhndrNSul+czGE8Y6SU90gpt0sptzc15ZaR/daZUZ45nJsAMhO7D/byzacO857z2rjtshUFO44Vsta7R0K05cGpDlp58+rKirLWRDxOe9L3kw/aat1MRBNlV0K/Qxciqy3kE+kCuqSUL+qvH0YTKr26KQr9sS9l/pKU7duBnjnG29OMM8sx8k6jt/Srhg4Ewnz58QOsaqrin997TkFr8mhZ6+VtzuoeDuW193WTr5K+ctVEhkIsqfPk9TvXrpvGys0vcrQ/QI3bQUNV7i2EzcScQkRKeQboFEIYMaJXAweAxwAjwmon8Kj+/DHgVj1K6yLAr5uingSuFULU6Q71a4En9ffGhBAX6VFZt07bV7pj5J1GbyUjwWjJ1uoJx+J86oevcto/wZdv3jyv9qSZUFflLOvGVOFYnB5/iGXz6AI5nZZqV1lrIvmomZVKm+6k7y5DIbKqqfQLLxpk2iThk8BDQggncAz4YzQB9BMhxG3AKeAWfe4TwI1ABxDU5yKlHBJCfBl4WZ/3JSnlkP7848D3ATfwC/0P4KszHCPvNHgno2da8lCFdCGZiMa5/Qd7ee7oIF+/5VwuWdVY8GPWeRwMB6MkEjIvGcpmo3MoiJSwvDF/Nv5mXyWvnBrO2/7MgpSSruEQF61syOt+k0JkJJjX/Rabjr5xrlo//0KoZiEjISKl3AdsT/PW1WnmSuCOGfZzH3BfmvE9wOY044PpjlEIkhmyY+GSEiLReII//69X+M3hfr723nN4X57DeWeirdZDPCHpHgnl1ZlqFk4MaBeu5fnURHyaJiKlLJu7UNAalAXCsbyG9wL4XA6qKyuSodblwPB4hIFAuGyc6qAy1pM06VVDzeoXiSckR/sDjE1M+iHGJqJ8/D9f4amDfXx5xyb+8PylC7aetS3aj+BIX3lmrZ8Y1LrO5VOINFdXEo4lGA3F8rZPM5Cs3luAm4mWGldZdYTMR4tqs5F7z88yo1VPKDNb97mJaJw7HnqF3W9NxhRcvb6ZjYt9PLy3i76xMF/asYkPX7x8Qde1pkXLWj/cG+CqPFUHNhMnBsepcTuoy6Pz08gV6R2boCaPUUzFxnB851sTAcqumdfr3X4ANrf5iryS/KE0EZ2W6koqbMJUkSChSJzP//S1KQIE4LcdA3zr6Q4W1bj4yccu5tYFFiBAMrrk5GB52asNTg4GWd6Q3zvrlmQ9KHNqu7nSOVRATcTnotdkN3bz4Y0eP+11bmo95RGZBUoTSVJht7Go1pXMvC02g4EwO+9/iTe6R/nsNWsB+JddWm+Lxz95GcsaqnDOs8fFfFmsJ4OVI8cHxtm6NL+l2pKaSBndWYNmzqpxO/LS8Gw6rTWV9I2FiSdkXpI+i82b3f55txUwG0qIpNBe6zGFJjI0HuED332BU0NB7t25nas3tCCl5MMXL+ON7tGkKanYLKpxJX0H5UQ4FqdnJMR7tuY3SMHIWu8rsx4ZnUP5q947nVafi1hCMhgIJ4VwqTI6EeXEYJBbti+Ze3IJocxZKbTXFb/g2+hElA/f+yInB4Pc95HzuXqD5m/Q6lU5uWxN4cN3M2VxrZvTZRQ5Y9A5FCIhybs5q6qyAm9lRdlpIh19AVYWKNrIiJQsB7/Im92aU31zmWkiSoik0F7noXdsgnCseI2D/vZnr3PozBjf/vC2Bcn3mA+LalyMhWNTIsbKgZNGZFYB+l83+yrLqlvfeDhG90iIdS2FESKtNboQKQO/yBuGU31x+TjVQQmRKbTXuZGSosWlv97l5/HXTvOJK1dz5bq0tSZNRbPP6INdXpnrJwbznyNi0FxdWVaayBG9DlShTKytZeRHeqPHz+IaFw156jBqFpQQScGw6xbDuR6LJ/jbR16n0evkT95euMKJ+aShyhAi5XNnDXBiYByfqyKvxQQNWnyusvKJHNa7W64tkBBp8FZit4myMGe93u1nU5mZskAJkSm01xev4NuPXu7k9W4/X7xpc0GiXAqBUSpmsNyEyOA4yxsLU9uoxaclz5VLefPDvWNUVthYWqCqBXaboNXnKvms9UA4xvGB8bKLzAIlRKbQ6nPpuSILq4nEE5Lv/t8xtiyp5cZzWhf02PPBKBVTfuas8bwWXkwlmbU+UR5Z64f7Aqxu9hY0/FYLeDFH6H2uHOgZRcrySjI0UEIkBbtNsLgILTmfOtjLycEgH7t8ZUnVVKqvMjSR8hEiE9E43cMhVhTAqQ6pHQ5L+87a4PCZsYKZsgza6zzJ9rulStKprjSR8qcYYb7/vaeTVp+LazeVjhYC4LDbqPU4ysoncqQ3QELC+tbCXBjLKVfEH4pyZnSi4EJkSb276FGT8+WNbj/N1ZU0V5d2rks6lBCZRnudm5OD4wtms+4bm+DZw/3ctGVxSWbktlSXV4G8t85osfyFEiJGtFE5ZPp39BlO9cJWpF1S50FK8/QVkVLyz784yB999wUe298z57VCSskLxwbZtiy/FRDMghIi0zinvZaBQIRTQwtjg/3P508SS0j+6IKFq8CbT1pryqtA3ltnxnA5bAXzibTVubHbRFnUHHtDT55bv6iwdn6jJlenSYTIo/t6+M6zx3jp+BCf+uGrPPTiqVnnHx8Yp8c/wSWrzZ33lStKiEzjohX1ALx0fGiOmfMnHIvz0IunuHp9S0ES2xaCRTUuTpdBIpjBoTNjrGmuLphW6LDbWFLn5vhA6ZeLeeXUMC2+ShbXFNZEU8zQ+3Q8vLeLZQ0eDn3lBi5YUc+XHj/Aa10jM85/5NVuhICr1ps/9ysXlBCZxsomLzbBgmgiv36rn8HxCB+6qDS1ENA0kYFAmEisNNsKT+etM2MFM2UZrGisKgsh8uqpEc5bUlfwYJAWnwuHXZjCud43OsFzRwfYsaUNu03wnQ9to97j5NM/2pc21N0fivKfL5zkynXNtNUWpr5YsVFCZBp2m6DRuzBZxY/u66bR6+SyElZzF9doWf4LZf4rJEbXuUI7ilc0ejk+sHB+t0IwEAhzaijI1mW1BT+W3SY0s6m/+ELk56+dJiHhpnMXA1BX5eRbf3Qe3SMhdt7/EiPBSPL/ur9zhB3/328ZnYjxqavXFHPZBUVV8U2DlhBW2OiZ0Ykou9/q44MXLqXCXrqy/OJVWl/tp9/qZXVzabf8PK7XzCpUeK/BikYPoWic3tFwsjZUqfHKSa1XfL7L5c9ElbOC8Ujxo7Me29fN5jbflO/6+cvruesD5/GxH+xly5d2AVoARd/YBItq3PzgtgvYsqTwwrZYlO7Vq4C0+AqvifzuyACRWIIbz1lU0OMUmiX1HlY1VfHS8eFiL2XenBgoXOHFVFY0ahegUjZpvXJqhAqbWLC8B4/TTqjIQuT4wDj7u/zsOLftrPcuX9OUfO6wCy5Z3cAnr1rDzz95mekLqc4XpYmkodnnYu/JYaSUBbP3PnOon2pXBeeVwR1Ke52nLMJ8j/WPYxNaXkIhWd6oRRsdHxhPanLF5vmjg3xj1yFuPq+ND164bM75r54aZtNiHy6HfQFWBx5nBcFIcbP8H92nOch/TzdlpeJyTN6Pf/MPt/Dut509p1xRmkgazm2vYTgY5fmjgwXZv5SSZw/38/Y1jSVtyjIwWx/sSCzBP//iIFf/yzM8uq874+2efquPt7XXUllR2Avj4ho3zgobxwcCBT1Opuw9OcwHv/cCL58Y5u8eeYNdB3pnnR8Ix9jXOcK2ZfULtEJwO+0Ei6iJSCl5dF8PF61oSGuCTL3ZLJXad/mi9K9gBWDHljZcDhu7Ds7+Y8qVQ71jnBmd4Iq15RHy16JHaEXj5ojQ+vLjB/jOs8foHQ3z2Z/szyhcu6MvwIHTo0mHaSGx2QTLGzwcHzBHMMK3nj5Co7eSPX//Ts5pq+EzP97Hmz3+s+ZF4wkGA2F+vr+HcCzB9ZsXrsKCp8hC5PVuP8cHxtmxZe7vh8+thIjlcTnsrG720tFXmDvFZw71A3D52qY5ZpYGrT4XUmKKZkuhSJwf7+nk/ecv4bkvXMWSeg+feGjvnE2NHn+tByHgXW9bGB+VFuZbfE1keDzCb48M8Ptb22j0VvLtD2+j2lXBrfe+lMxI94c0rfzif36abV95ii/87HXOW1rL9gXMwPY47ZwaCrL3ZHF8bz97pRuHXXDD5rm/HzVKiCgA1jRXF0yIPHuon/Wt1SUbmTMdIxnshAkcxc8f0wIW3vW2RfhcDr7z4W0EI3E+cv9LyZIm05FS8vP9PVy4oj7ZjrXQrGj0cmooWHTt7ZdvniGWkPyebsNvq3Xzn39yIUIIbvn287znP37HuV/8FR/47gsMBMJcvraJnRcv496d52NbwDI9bofmvn3v3c8t2DENukdCPLy3ixs2L6Imgx4zPpe1XM3WOtssWN3s5ZFXuxmbiFKdRxtnIBxjz8khPnpZaTSeygSjR8L+Ln/RSzs8c6gft8POBXrlgbUt1dz9oW385Y/38e67fsuyBg91HifNPq0Y3js3tFDpsHG0f3xB/yermqqIxiWnhoKsKlB/8kzYfbCPZQ0eNqW0bF3Z5OWhP7mQP/+vVxgcj/Cpq9ewtF6bs66lekGFh4HHuTAO/OlEYgn+9IE9CODPr1qd0Tb5vF6UAkqIzIARB360fzyvMd7PHOojGpcl0f42U+qqnCxr8LC/c+bSDwuBlJKn3+rj0tUNU5zj71jbxFOfeQfffvYo3cMhhsYjHDozxrOH+vn+cydw2m00ep2857z2BVtr8vvVFyiqEDnUO8qWNFnn61qr2fWZdxRpVWeTKrcC4RjeyoW5dD3w3AkOnB7lng9vyzgJ1VlhLQOPEiIzYPzIO/oCeRUij+8/TYuvkvOXL1xky0KwodXH4d6xoq7haP84XcMh/uwdq856r77Kyd/euGHK2EQ0zlf+9wBH+8b53A3rcS/g3e4q4/vVH+DaBTvqVMbDMTqHQvzBtiVFWkHmBMKTTvXj/eOc0174/JTh8Qh3PX2EK9Y1ZdSm4eE/u5h9Rb6RKgZKiMzAsnoPDrvgSF9+L4ydw0E2La4pybLvs7GmxcuvDpwhHIsXPER2Jp451AfAFesyC1hwOex85eZzCrmkGfG5HLT4Kgvmd8uEV05pTup1Ba4Vlg8C4Wjy+bGBwIIIkZ++0sXYRIzPXb8+o/nbl9ezvcxuDjPBWnpXFlTYbaxs9NLRm98feTASp2qBVPGFZHWzl4SEE0UMW33mUD9rmr201xWm33e+Wd3s5WgRhcgDz52k0VtZElGCa5onBd1CZPpLKfnRy52ct7SWDQUudV/qKCEyC6ubvXT05/dHPh6O4a0szp16ITHMf/nW3DJlPBzjpeNDXFlC5bZXNXk52l+cQozhWJzfdQzwrnNaFyzrfD7cdtkKHv/kZbQvUBn9l08M09EX4P3nm9/UV2wyEiJCiBNCiNeFEPuEEHv0sXohxC4hxBH9sU4fF0KIu4QQHUKI14QQW1P2s1Off0QIsTNlfJu+/w59WzHbMRaK1c1aGGY+a/YEI3E8zvLTRFY1eRFCay9bDJ47OkgknuCKErirNljd7CUQjhW82Gc6Xj01Qiga57I1pfF52fQ6XSubCpe/lcp3nj1KfZUzbYkTxVSy0USulFJukVJu119/HtgtpVwD7NZfA9wArNH/bgfuBk0gAHcCFwIXAHemCIW79bnGdtfPcYwFYeNiH1LCgdNnZ+/mwjd2HSYQjlFVpHDFQuJy2Fla78m75pYpvz7UR5XTXlI26dVNk8EbC83vOgaw2wQXriydzwtgTbOXo/0BEonCaW9n/BP8+pBWYbscb/jyzXzMWTuAB/TnDwA3p4w/KDVeAGqFEIuA64BdUsohKeUwsAu4Xn/PJ6V8Xmp6/YPT9pXuGAvCeUu1qKxXT80/4iIaT3DX7iMAeMrQJwLaDzzfPqRMkFLy7KF+LlvTWFLhlZMRgAtvAvzN4X7Oba8puTpPa1u8TEQTdBawy+HPXu0iIeF92xYu5LuUyfQXJ4FfCSH2CiFu18dapJSnAfRHwxjdBnSmbNulj8023pVmfLZjTEEIcbsQYo8QYk9/f3+GpzQ3zdUu2uvcySiW+fBmz2S2dDlqIgCrm6s5NhAgtsBZ2L2jYbpHQly80hwVcTOlqboSn6uCwwusiZwaDLK/y891GYStmo2VuvZ2okA96qWUPLy3iwuW17OsoTRbVi80mQqRS6WUW9FMVXcIIS6fZW662FWZw3jGSCnvkVJul1Jub2rKr433vKV1edFEUmv+lKuKvLrZm8zCXkhePqEVWCy1KBohBBsW+Th4On05lkLxvd8eo8ImStLe79LDx8PRwhRjfLVzhGP940oLyYKMhIiUskd/7AMeQfNp9OqmKPTHPn16F5Aa0tAO9Mwx3p5mnFmOsWBsXVrLaf8Ep+fZmnPvyclKssUq4VBo1ujmmYVMOjzWH+CTP3wVgPWtpSVEQBN8b50eI15AG38qnUNBfvjSKW7ZvoTFJdjzu8Ku3XPGCvR5Pby3C7fDzo0LVIizHJhTiAghqoQQ1cZz4FrgDeAxwIiw2gk8qj9/DLhVj9K6CPDrpqgngWuFEHW6Q/1a4En9vTEhxEV6VNat0/aV7hgLxtvaNb/IG9253y1KKdlzYlITKdQPoNisbanGJuDg6YUTIt98SvMzvettmRXHMxsbF/sIReOcHFyY4pV37T6CEIJPXZ1ZHSiz4dD77xSicOVENM7P9/dww+bWBSurUg5k8km1AI/oUbcVwH9JKX8phHgZ+IkQ4jbgFHCLPv8J4EagAwgCfwwgpRwSQnwZeFmf9yUppXF7/nHg+4Ab+IX+B/DVGY6xYGxYpF0YX+/2c83Glpz20TUcoi+lTHqoQKp4sXE77axorOLAAplnDvSM8vP9Pdxx5Sr++rrMsorNxkbdBHfw9FjS3l8ojvUH+OkrXXzkkhUsqik9LQS01rMAsXj+b8R+daCXsYmYMmVlyZxCREp5DDg3zfggcHWacQncMcO+7gPuSzO+B9ic6TEWEo+zglVNXt7szj3M1/CHfPtD27jvt8e5bmPpOTQzZePiGl5ZoJ4P39h1CJ+rgtvffnatrFJhdbOXCpvgwGl/wXuZfPOpI7gcdj5xZel+XkYn0Fgi/5rIw3u7aKt1c1GJBWgUm9KJhywi57TV8EaaTm+ZsufkEFVOO9dsbOEnf3ZxSZpdMmXjIh/dIyH8oejck+fBK6eGeepgHx97x6qS/jxdDjurmrwFNwEePK1pbX986XIavZUFPVYhceg156J51kQ6h4L835F+3ru1rSil7ksZJUQyYFNbDb2jYfrGcusjvvfkCOctrSu7oovp2LBIq3FU6Iijrz95iEavk49csrygx1kINi72caCnsJ/Xv/zqMNUlrrXBpCbS0RfIa7mYh148hQDef8HSvO3TKighkgFG06XXu7LXRsYmohw6M8q2BWwlWkw26s2NCnlR/F3HAM8dHeQTV6wui2KWGxZVc2Z0gqHxSEH2/+KxQZ462MvHLl9Z0lobTPpEvv/cCb71dEde9jkejvGTPZ1cs7GlJCPWio0SIhlwTptWuj2XfJFXT42QkLB9uTWESHO1i0ZvZcE0kXhC8tVfvMXiGhd/dGF53DWe225URsi/LymekHzp8QMsrnFx22Ur877/hcaIzgL4zxdO5mWf3372KEPjET6Wpg+NYm6UEMkAt9POxkU+Xjw+mPW2e08OYxPktbGV2dmwqHpKhn4++a8XT/J6t5/P3bC+JKrPZsK5S2px2AUvn8i/EPnp3i7e7Bld8KZbhaIixSScGvGYK72jE9zzm2Ps2LKYrUutcaOXb5QQyZCr1jez5+Rw1n6RvSeHWdfqs1Tf5XPaajjcO5bX6scG9z93gu3L6ripBLOtZ8LlsLNpcc2UhNR8MBGN8y+7DnHe0tqy+bym+xUD4di89nfPb44RS0g+e826ee3HyighkiE3nrMIKeHJN85kvE0snuDVU8Nst4g/xGDbsjpiCcn+rvy2Co0nJL3+Cc5przmrJ3ipc/7yOvZ3+QnH8id47/3tcXpHw3z++vVl83lNP4/982hH2zkU5AfPn+Q957WxtKE0GpmZESVEMmRti5dVTVU88XrmQuRQ7xjjkbhl/CEGhllgbx7zRUKROGv//heMR+I0VZduiOpMnL+8nkgskZc6bQBPvnmGb+46zPWbWrmwjPMeftsxkPO2X//VIWw2+My1a/O4IuuhhEiGCCG48ZxFvHh8kIFAZrZY4yJqNVtrXZWTVU1VeRUi//TEwWR9qaYSznOYiYtXNVBhEzx7eH5VqPtGJ/jcw6/xsR/sZdNiH19739vytELz0VDlZPfB3py23d85wqP7erjtstLN3jcLSohkwY3nLCIhtbu8TNhzYpgWXyXtddb7km5fVs/ek8N5KSz40Isn+UFKJE5jGWoi1S4H25bV8cyh3IXIC8cGede3fssjr3az8+Jl/PhjF1PjLl9f3LvetojDvQFOZVkWPpGQ3PnYmzRVV/JnKiJr3ighkgXrW6tZ0VjFLzI0ae09Ocz2ZfVlY4/OhsvWNOIPRecdtvpcxwB3PvomV65rSpZ6L9fieFesa+bg6VF6R7ML3nij28/b/9/TvP+eF/BWVvDzT17GF3dsLpvotZm48RytTMxTWWoju9/qY1/nCH9z3TpLBbwUCiVEskAzabXy/LHBORPDOoeCdI+ELJNkOJ13rGuiwiZ46mDu1fuPD4zz8YdeYUVjFXd94Dy+eNMmljV4WN9anceVmocr1mm9cLIxaU1E43z8ob1EY5J/ePdGfv7Jy1hXpp/PdNY0e1nb4s1aiPx0bxeNXie/f17b3JMVc6KESJbcsHkR8YTkV3OYtIwv9lXr0zZjLHt8LgcXrKjP2WbdPRLiw/e+iE3AvTvPp1rf37N/fWXZ3j2ub62m1efimUOZC94fvnSKzqEQ3/iDc/noZSvKVktLh7PCxjvWNvHS8aGMzaYjwQhPv9XHTee2JUuoKOaH+hSzZNNiH0vrPfzv66dnnfc/r3azrqWa5Y3WbbF59YYWjvRlb7P+9aE+rv/mbxgJRnngoxdYJvxSCMGV65t59lA/Exm0C5iIxvmPZ45y0cp6LlnduAArNBfOChuLa93EEpLRDAt+/upAL5F4QmkheUQJkSwRQnDzlsX8tmOAjr70lVcPnRljf5efPzh/Sdr3rcI7N2ha2K4stBF/KMpnf7Kftjo3/3PHpcmmYFbh+s2tjEfiPP3W3NrIQy+eon8szF++05ohqk67jfoqJwBDwczqjj1zqI9Wn4vNbaXXBdOsKCGSAzsvWU6Vs4LP/fR1Ymk6rD3yajd2m2DHlvLIEs6VZQ1VbFzk49F93Rlv8x/PdDAcjPD1W85ldXNhmzSZkUtXNdBe5+b+3x2fdV4oEufuZ45yyaqGss4DmQ0hBHUeTYgMZ1C8MhZP8H9HBnjH2iZLBrsUCiVEcqDBW8k//v5m9p4c5l/19qwGiYTksX3dXL6msaT7NuSLW7a381qXP6OCjMPjER587iQ3b2ljs1452WpU2G185JLlvHximNdmyfh/8PkTDATCfOYaa2ohBoYQyaQC8iunRhibiCUDGBT5QQmRHNmxpY1btrXz7890TCkR/+tDffT4J7hZ2VwBuHlLG067jR+/3Dnn3AeeP0EoGucTV1g7dv8Pz1+Ct7KC7/5fem1kdCLK3c8e5Yp1TWxfXr/AqzMXdVVakMVwBuasZw71UWETXLrGev6jQmKdUI4C8Pfv3sgzh/u57YGXOXdJLcFIjJdPDLO03pOMYbc6dVVO3vW2Rfz3nk7+8pq1Mya/+YNR7v/dCd65oYU1LdYIUZ2JapeDD120jG8/e5T3bm3jinVTI/y++5tjjASj/NW1qmhg0icyPrdj/ZlD/WxdVoevTKP7ioXSROZBjdvB/R85nzUtXt46M0ogHOfSVQ18+0PbpvQ9sDp/8vYVjEfis/Z/+PdnOhidiPJX11nbPGPw6XeuYfXrdeEAAAh0SURBVF1LNZ/+8T5eOq5V940nJD966RTffvYoN5272LImv1TcDjveygpO+0OzzuscCnLg9KhlQ+4LidJE5snmthoe+pOLir0MU7NpcQ1XrW/mW08f4frNraxqmuow7xwK8v3fneB9W9tZ36qiZkArD3/Prdv44/tf5oPfe4H3bm3nxOA4Lxwb4tz2Gr5406ZiL9EUCCFY1eyloy8w67xf6tW3b9ysLAT5Rt0uKxaEf37PObgcdj79o31EYlMj2r7+q0MIoaqpTmdZQxWPfOJSrl7fwmP7e+gaDvHX163jkU9cSp1uxrEiT376cn72iUuSr9c2ezkyhxB54o3TbG7zWSbnaCFRmohiQWjxufjyjs188oev8rVfvsVndYFx1+4OHt3XwyevWq2qqaahxuPg2x/eVuxlmIrpZV3WtVbz33u76BwKsqT+bCHRMxLi1VMj/PV1yodUCJQQUSwY121qBbRmST99pQuAkWCUP9y+hL+4ek0xl6YoYa7b1MpX/vcgj+3v4Y4rV5/1/hN6dYkbNrcu9NIsgRIiigXDWWHj3W9bxJs9o2xuq8FhE3zwoqVsW2btMFXF/FhS7+GSVQ08+PwJPnrpiim95GPxBD944SRbl9ayssl6yasLgRIiigXlWx84j3hCquJ3irzy6Xeu5Q++8zz/9MRBvnzz5uT4j/d0cnIwyBdu2FDE1ZU36pesWFCEEEqAKPLOBSvq+dO3r+AHL5zkrt1HiMUTPNcxwD/+70EuWlnPdZtair3EskVpIgqFoiz4m+vXMxiI8I1dh7nvd8fxh6Ksafbyb+8/T9XKKiBKiCgUirLAYbfxL39wLtdvbuXJN3tZVOPiTy9fWdYtgs2AEiIKhaJsEEJw7aZWrt2kIrEWCmWcVigUCkXOZCxEhBB2IcSrQojH9dcrhBAvCiGOCCF+LIRw6uOV+usO/f3lKfv4gj5+SAhxXcr49fpYhxDi8ynjaY+hUCgUCnOQjSbyF8DBlNdfA74ppVwDDAO36eO3AcNSytXAN/V5CCE2Au8HNgHXA/+hCyY78O/ADcBG4AP63NmOoVAoFAoTkJEQEUK0A+8Cvqe/FsBVwMP6lAeAm/XnO/TX6O9frc/fAfxIShmWUh4HOoAL9L8OKeUxKWUE+BGwY45jKBQKhcIEZKqJ/CvwN4BROa8BGJFSxvTXXYDRhakN6ATQ3/fr85Pj07aZaXy2YygUCoXCBMwpRIQQ7wb6pJR7U4fTTJVzvJev8XRrvF0IsUcIsae/vz/dFIVCoVAUgEw0kUuBm4QQJ9BMTVehaSa1QggjRLgd6NGfdwFLAPT3a4Ch1PFp28w0PjDLMaYgpbxHSrldSrm9qUn1T1YoFIqFYk4hIqX8gpSyXUq5HM0x/rSU8oPAr4H36dN2Ao/qzx/TX6O//7SUUurj79ejt1YAa4CXgJeBNXokllM/xmP6NjMdQ6FQKBQmQGjX6gwnC3EF8FdSyncLIVaiaSb1wKvAh6SUYSGEC/gBcB6aBvJ+KeUxffu/Az4KxIBPSyl/oY/fiKbd2IH7pJT/qI+nPcYca+wHZu7DOjuNaBpQOaHOqTRQ51QalPM5LZNSZm3KyUqIlDtCiD1Syu3FXkc+UedUGqhzKg3UOZ2NylhXKBQKRc4oIaJQKBSKnFFCZCr3FHsBBUCdU2mgzqk0UOc0DeUTUSgUCkXOKE1EoVAoFDmjhIjOTJWEzY4Q4j4hRJ8Q4o2UsXohxC69+vEuIUSdPi6EEHfp5/iaEGJr8VaeHiHEEiHEr4UQB4UQbwoh/kIfL+VzcgkhXhJC7NfP6Yv6eNaVsM1GPqp7mwkhxAkhxOtCiH1CiD36WMl+9wCEELVCiIeFEG/pv6uL83lOSoig/RCYuZKw2fk+WlXkVD4P7NarH+/WX4N2fmv0v9uBuxdojdkQAz4rpdwAXATcof8vSvmcwsBVUspzgS3A9UKIi8iyErZJmVd1b5NypZRyS0rYayl/9wD+DfillHI9cC7a/yt/5ySltPwfcDHwZMrrLwBfKPa6slj/cuCNlNeHgEX680XAIf35d4APpJtn1j+0KgXXlMs5AR7gFeBCtASvCn08+R0EngQu1p9X6PNEsdee5lza9QvQVcDjaPXuSv2cTgCN08ZK9rsH+IDj0z/rfJ6T0kQ0ZqokXKq0SClPA+iPzfp4SZ2nbvI4D3iREj8n3eyzD+gDdgFHyb4SttnIR3VvsyGBXwkh9gohbtfHSvm7txLoB+7XzY7fE0JUkcdzUkJEI+OKwSVOyZynEMIL/BStPM7obFPTjJnunKSUcSnlFrS79wuADemm6Y+mPyeRv+reZuNSKeVWNLPOHUKIy2eZWwrnVAFsBe6WUp4HjDNpukpH1uekhIjGTJWES5VeIcQiAP2xTx8vifMUQjjQBMhDUsqf6cMlfU4GUsoR4Bk0f0+2lbDNRL6qe5sKKWWP/tgHPIIm8Ev5u9cFdEkpX9RfP4wmVPJ2TkqIaKStJFzkNc2H1ErK0yss36pHYFwE+A2V1iwIIQRwL3BQSvmNlLdK+ZyahBC1+nM38E4052a2lbBNg8xfdW/TIISoEkJUG8+Ba4E3KOHvnpTyDNAphFinD10NHCCf51Rsx49Z/oAbgcNotuq/K/Z6slj3D4HTQBTtLuI2NFvzbuCI/livzxVoUWhHgdeB7cVef5rzuQxNfX4N2Kf/3Vji5/Q2tCrUr6FdlP5BH1+J1g6hA/hvoFIfd+mvO/T3Vxb7HOY4vyuAx0v9nPS179f/3jSuA6X83dPXuQXYo3///geoy+c5qYx1hUKhUOSMMmcpFAqFImeUEFEoFApFzigholAoFIqcUUJEoVAoFDmjhIhCoVAockYJEYVCoVDkjBIiCoVCocgZJUQUCoVCkTP/P69DzPewfFdLAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "data.plot(figsize=(6,4))\n",
    "# print(data)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAs4AAAHiCAYAAAD4akr8AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzt3X+43VddJ/r3JydNm1JioE0qoekPIXZaeDRgBsQfM9Hi3NYflHFQqYrVW60+ijojV0XkQcWRcWauP4YZ7txhKAPWEQYYf3S0XHSKVccr3IYS1LbGhtqSkNqU0hBKQw85Z90/zg6cnJ6TrHP2ztnnJK/X8+Q5e3/3+u619srae7/32uv73dVaCwAAcGJrxt0AAABYDQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzQIeqek1VvaWz7Nuq6l+e6jatdFX1vVX1v4bY/71Vdf0o2wQwDMEZOC1U1f1VdaSqHquqh6rqv1TVeUu8r51VtX/2ttbaG1pr3z+a1n6+jlZVP7XI/X6+qn5zVO1YKeZ7XK21a1prbx9XmwDmEpyB08m3tNbOS/L8JP8wyWsXewdVtXbkrZrf9Uk+Ofi7otWMNSfbBnC686IHnHZaax9P8t4kz02Sqvq+qrqnqj5dVfdV1Q8eK3tsdrmqfrqq/j7JOwb7bhnMXj9WVVvmzohW1bur6u+r6lNV9adV9Zze9lXVuUleluRHkmyrqh1z2zOn/P1V9eKqujrJa5J8x6BdHxncvqWqbqmqT1bV3qr6gVn7TgyWmXx08Pg/VFVbB7d9VVXdMXgMd1TVV83a7/aq+qWq+vMkjyf5kgW2fVFV3VRVD1bVx6vqX1bVxAKP+99V1b6qOjxox9cOti/0uG6vqu8fXF5TVa+tqgeq6mBV/UZVfdHgtksHs/fXV9XHquoTVfWzvf8fAL0EZ+C0MwiG35jkw4NNB5N8c5INSb4vya9V1fNn7fLFSZ6e5JIk35PkmiQHWmvnDf4dmKea9ybZlmRzkjuT/NdFNPGfJXksybuTvG9Q50m11v6fJG9I8t8G7frywU3vSLI/yZbMBPI3VNVVg9t+Isl1memPDUn+9ySPV9XTk/xBkjcmOT/Jryb5g6o6f1aVr0hyY5KnJnlggW1vT3I0ybOTPC/JP0my0JKWO5Jsz0xf/1aSd1fVOSd4XLN97+Df1yX5kiTnJfkPc8p8TZLLk1yV5HVVdcUC7QBYEsEZOJ38blUdSvK/kvxJZsJYWmt/0Fr7aJvxJ0n+MMnXztpvOsnPtdaeaK0d6amotfbW1tqnW2tPJPn5JF9+bAa0w/WZCYlTmQmQ11XVWZ37HmfwIeFrkvx0a+2zrbXdSd6SmYCbzITY17bW9gwe/0daa48k+aYk97bWbm6tHW2tvSPJ3yT5lll3/7bW2l2D2z83d1tmAvA1Sf55a+0zrbWDSX4tycvna2tr7Tdba48M7u9XkpydmaDb47uS/Gpr7b7W2mNJfibJy+csrfmF1tqR1tpHknwkyXwBHGDJBGfgdPLS1trG1tolrbUfPhaCq+qaqvrAYCnDoczMvl4wa7+HW2uf7a1ksPzhlwfLHw4nuX9w0wUn2O3YvlszM2t6bIb695Kck5kguxRbknyytfbpWdseSPLMweWtST66wH4PzNk2e78k2TfPfrO3XZLkrCQPVtWhQd/+p8zMwj9JVb1qsGTmU4OyX5SOPlugvQ8kWZvkwlnb/n7W5cczMysNMDKCM3Baq6qzk/z3JP9nkgtbaxuT3JqkZhVrc3abe32u70xybZIXZyb8XXqsuo4mvSIzr73/Y7Cm+r7MBOdjyzU+k+TcWe2fSLLpBG07kOTpVfXUWdsuTvLxweV9SZ41TzsOZCb4zjZ7v/nqmrttX5Inklww+MCysbW2obX2pPXeg/XMP53k25M8bfD/8Kl8oc9O1udz23txZpaIPHSS/QBGRnAGTnfrMrMk4OEkR6vqmsyswz2Rh5Kcf4KlF0/NTGB8JDMh9w2LaM/3JPmFzKz1PfbvnyX5psH64r9Nck5VfdNg+cZrB+2f3bZLj53RorW2L8n/m+RfVdU5VfVlSW7IF2a035LkF6tq2+BMGF82qOfWJF9aVd9ZVWur6juSXJnk93sfSGvtwcwse/mVqtowOIDvWVX1j+cp/tTMBN2Hk6ytqtdlZs31vI9rHu9I8i+q6rKaOc3gsTXRR3vbCzAswRk4rQ2WMPxYkncleTQzs8W3nGSfv8lMULtvsARhy5wiv5GZpQIfT3J3kg/0tKWqvjIzs9Nvaq39/ax/tyTZm+S61tqnkvxwZgLvxzMzAz37LBvvHvx9pKruHFy+bnC/B5L8TmbWa//R4LZfHTz2P0xyOMlNSdYP1jl/c5JXZeYDwE8l+ebW2id6Hsss35OZDyd3Z6Z/35PkGfOUe19mDqj828z03Wdz/LKP+R7XbG9NcnOSP03yd4P9f3SRbQUYSrV2sm/HAAAAM84AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAd1p68yHhccMEF7dJLLx13MwAAOM196EMf+kRrbdPJyq3Y4HzppZdm165d424GAACnuap6oKecpRoAANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoMJLgXFVvraqDVfXXC9xeVfXGqtpbVX9ZVc8fRb0AALBcRjXj/LYkV5/g9muSbBv8uzHJfxxRvSM1Nd1y2z0P5Y233Zvb7nkoU9Nt3E0CAGCFWDuKO2mt/WlVXXqCItcm+Y3WWkvygaraWFXPaK09OIr6R2FquuUVN30wu/cdypHJqaxfN5HtWzfm5htemIk1Ne7mAQAwZsu1xvmZSfbNur5/sG3FuH3PwezedyiPT06lJXl8ciq79x3K7XsOjrtpAACsAMsVnOebsn3SOoiqurGqdlXVrocffngZmvUFdx04nCOTU8dtOzI5lbsPHF7WdgAAsDItV3Den2TrrOsXJTkwt1Br7c2ttR2ttR2bNm1apqbNeM6WDVm/buK4bevXTeTKLRuWtR0AAKxMyxWcb0nyPYOza3xlkk+tpPXNSbLz8s3ZvnVjamoyadM5d7DGeeflm8fdNAAAVoCRHBxYVe9IsjPJBVW1P8nPJTkrSVpr/3eSW5N8Y5K9SR5P8n2jqHeUJtZUbr7hhXnRt96Qyadszq+89l9k5+WbHRgIAECS0Z1V47qT3N6S/Mgo6jqVJtZUzj10X849dF+uuuLCcTcHAIAVxC8HAgBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6DCS4FxVV1fVnqraW1Wvnuf2i6vqj6vqw1X1l1X1jaOoFwAAlsvQwbmqJpK8Kck1Sa5Mcl1VXTmn2GuTvKu19rwkL0/yfw1bLwAALKdRzDi/IMne1tp9rbXJJO9Mcu2cMi3JhsHlL0pyYAT1AgDAshlFcH5mkn2zru8fbJvt55N8d1XtT3Jrkh+d746q6saq2lVVux5++OERNA0AAEZjFMG55tnW5ly/LsnbWmsXJfnGJDdX1ZPqbq29ubW2o7W2Y9OmTSNoGgAAjMYogvP+JFtnXb8oT16KcUOSdyVJa+0vkpyT5IIR1A0AAMtiFMH5jiTbquqyqlqXmYP/bplT5mNJrkqSqroiM8HZWgwAAFaNoYNza+1oklcmeV+SezJz9oy7qur1VfWSQbFXJfmBqvpIknck+d7W2tzlHAAAsGKtHcWdtNZuzcxBf7O3vW7W5buTfPUo6gIAgHHwy4EAANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAh5EE56q6uqr2VNXeqnr1AmW+varurqq7quq3RlEvAAAsl7XD3kFVTSR5U5JvSLI/yR1VdUtr7e5ZZbYl+ZkkX91ae7SqNg9bLwAALKdRzDi/IMne1tp9rbXJJO9Mcu2cMj+Q5E2ttUeTpLV2cAT1AgDAshlFcH5mkn2zru8fbJvtS5N8aVX9eVV9oKquHkG9AACwbIZeqpGk5tnW5qlnW5KdSS5K8mdV9dzW2qHj7qjqxiQ3JsnFF188gqYBAMBojGLGeX+SrbOuX5TkwDxlfq+19rnW2t8l2ZOZIH2c1tqbW2s7Wms7Nm3aNIKmAQDAaIwiON+RZFtVXVZV65K8PMktc8r8bpKvS5KquiAzSzfuG0HdAACwLIYOzq21o0lemeR9Se5J8q7W2l1V9fqqesmg2PuSPFJVdyf54yQ/2Vp7ZNi6AQBguYxijXNaa7cmuXXOttfNutyS/MTgHwAArDp+ORAAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAECHkQTnqrq6qvZU1d6qevUJyr2sqlpV7RhFvQAAsFyGDs5VNZHkTUmuSXJlkuuq6sp5yj01yY8l+eCwdQIAwHIbxYzzC5Lsba3d11qbTPLOJNfOU+4Xk/ybJJ8dQZ0AALCsRhGcn5lk36zr+wfbPq+qnpdka2vt9090R1V1Y1XtqqpdDz/88AiaBgAAozGK4FzzbGufv7FqTZJfS/Kqk91Ra+3NrbUdrbUdmzZtGkHTAABgNEYRnPcn2Trr+kVJDsy6/tQkz01ye1Xdn+Qrk9ziAEEAAFaTUQTnO5Jsq6rLqmpdkpcnueXYja21T7XWLmitXdpauzTJB5K8pLW2awR1AwDAshg6OLfWjiZ5ZZL3Jbknybtaa3dV1eur6iXD3j8AAKwEa0dxJ621W5PcOmfb6xYou3MUdQIAwHLyy4EAANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAh5EE56q6uqr2VNXeqnr1PLf/RFXdXVV/WVW3VdUlo6gXAACWy9DBuaomkrwpyTVJrkxyXVVdOafYh5PsaK19WZL3JPk3w9YLAADLaRQzzi9Isre1dl9rbTLJO5NcO7tAa+2PW2uPD65+IMlFI6gXAACWzSiC8zOT7Jt1ff9g20JuSPLeEdQLAADLZu0I7qPm2dbmLVj13Ul2JPnHC9x+Y5Ibk+Tiiy8eQdMAAGA0RjHjvD/J1lnXL0pyYG6hqnpxkp9N8pLW2hPz3VFr7c2ttR2ttR2bNm0aQdMAAGA0RhGc70iyraouq6p1SV6e5JbZBarqeUn+U2ZC88ER1AkAAMtq6ODcWjua5JVJ3pfkniTvaq3dVVWvr6qXDIr92yTnJXl3Ve2uqlsWuDsAAFiRRrHGOa21W5PcOmfb62ZdfvEo6uH0MTXdcvueg7nrwOE8Z8uG7Lx8cybWzLdcHgBgZRhJcIbFmJpuecVNH8zufYdyZHIq69dNZPvWjbn5hhcKzwDAiuUnt1l2t+85mN37DuXxyam0JI9PTmX3vkO5fY/l7wDAyiU4s+zuOnA4Ryanjtt2ZHIqdx84PKYWAQCcnODMsnvOlg1Zv27iuG3r103kyi0bxtQiAICTE5xZdjsv35ztWzempiaTNp1zB2ucd16+edxNAwBYkODMsptYU7n5hhdm073/Ixv3/3n+/XXPc2AgALDiOasGYzGxpnLuofty7qH7ctUVF467OQAAJyU4s2TOxQwAnEkEZ5bEuZgBgDONNc4siXMxAwBnGsGZJXEuZgDgTCM4syTOxQwAnGkEZ5bEuZgBgDON4MySOBczAHCmcVYNlsy5mAGAM4kZZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA7OqsGqMzXdcvueg7nrwOE8Z8uG7Lx8s9PgAQCnnODMqjI13fKKmz6Y3fsO5cjkVNYPfnjFOaQBgFPNUg1Wldv3HMzufYfy+ORUWpLHJ6eye9+h3L7n4LibBgCc5gRnVpW7DhzOkcmp47YdmZzK3QcOj6lFAMCZQnBmVXnOlg1Zv27iuG3r103kyi0bxtQiAOBMYY3ziP3FRx8ZdxOW1eEjn0uytMe9lH3PWTuRyy54Su762CeSibU5+6y1ueyCp+SctRNnXN8DwOnmRc86f9xNOCEzzqwqa9ZUXnPNFTnv7t/N+r/7s/zY12/La665ImscGAgAnGKCM6vOmjWVdY/szfoH/jzPv+RpQvMpNj3dcucDj+a379yfOx94NNPTbdxNAoCxsFQDWND0dMsb3ntP9h58LJNHp7Nu7Zo8e/N5ZvkBOCOZcQYWtHvfoew9+FieODqdluSJo9PZe/Cx7N53aNxNA4BlJzgDC7r/kc9k8uj0cdsmj07n/kc+M6YWAcD4CM7Agi49/ylZt/b4l4l1a9fk0vOfMqYWAcD4CM7AgrZv3Zhnbz4vOTqZtOmcPVjjvH3rxnE3DQCWneAMLMjp/wDgCwRn4ISc/g8AZjgd3Wlgerpl975Duf+Rz+TS85+S7Vs3CjcAACMmOK9yzrMLQI8zcZLlTHzMnFqC8yo3+zy7yfHn2X3+JU8bc+sAWAmGnWRZjQF03BNLq7HPODnBeZU70Xl2Bef5eTEDzjTDTLKMO4Au1TgnllZrn3FygvMqd+w8u0/MCs/Os7uwM3HWZbXS1zA6w0yyrNZvNsc5sbRa+4yTE5xXuWPn2b3rY59IJtbm7LPWOs/uCZyJsy6rkb6G0RpmkmW1frM5zoml1dpnnJzT0a1yzrO7OMP8hPTs0N1yfOhmtPQ1jNYwP2a0Wn9BdJw/4LRa+4yTE5xPA86z22+YF7NhQjeLo69htIaZZFmtvyA6zoml1dpnnJzgzBnlTJx1WY3O1L6enm6584FH89t37s+dDzya6ek27iZxGlnqJMu4v9kc5nkxromlcfcZp441zpxRjr2Y/eCPvypT512YV/7Qjd0HnVlPvnxWc18v9aBG67pZyY4F0DyyN8+/5KeXrd7V/LwYV59xagnOnHGW+mI2TOhmcVZrXw/zJu8ofHiycZ9Szpl9mEtwhkUwg7B8VmNfD/Mm7yh8eLJxPS9W80w3p5Y1zgAjMsxBjWfqum4WZ7Wug19qu8f1vHBmHxYykhnnqro6yb9LMpHkLa21X55z+9lJfiPJVyR5JMl3tNbuH0XdACvFMOeNHXZdt6+VF2ec/XWmrYMfpt3jOt7BN0AsZOjgXFUTSd6U5BuS7E9yR1Xd0lq7e1axG5I82lp7dlW9PMm/TvIdw9bN8LzZwugM8yY/zLru1RqoxmWc/XUmroMfpt3jOt7Br/KykFEs1XhBkr2ttftaa5NJ3pnk2jllrk3y9sHl9yS5qqq8mo/ZsRfwN77/3rznQ/vzxvffmze8955V89UfrDTDnoJqqafO8rXy4oyzv4ape7We33zYdo/jlHLjPg/zal2Scyao1ob7z6iqlyW5urX2/YPrr0jywtbaK2eV+etBmf2D6x8dlPnEQvf79EuuaN/wmrcO1bal2P2R3UmS7V++fUn7H/7s50bZnG733v3XSZJtVz63e59Pf/ZoPn7oSGYPgarkmRvX56nn9H0ZsZR6x73vuOpureWxJ6by2c9N5ZyzJnLe2RNZLZ8fh+3v1VbvsJZ7fD386Sfyiccmn7R903nrcsFTz150GxZrtY3tcfbXMHWv1tfs1dru1lr+du99ycS6bNnyjCWN66XU3VrLxz55JEc+N5XWZvpq/VkTufjp61f082pYx15H0pJzz57IxvVnLevjfdcPfdWHWms7TlZuFGuc53tUc9N4T5lU1Y1JbkyS857xrOFbtgRLDczJTOiemmpjCYJL2e+zgyflbK0lT3xuqvvFbJhAM659x1H3sRfCx5/4XJJKralFvxCO88PCMP21Gusddv/lHl/nnDWRqjwpmJx91kTX/sMEhHGP7aXsO2x/javu886eyPqzJp4UqM47u7/d43jdXa3trqpcvm24LLKUuh97YurzfZXMjJUjn5vKY0/0vzevtkmtua8ja9ZUzjt7bf7BFz91xX1YGMWM84uS/Hxr7X8bXP+ZJGmt/atZZd43KPMXVbU2yd8n2dROUPmOHTvarl27hmrbctu5c2cOH/lc3vRbtyx63x/5zpckyZL2Xao7H3g0b3z/vcet4Tp77Zr82NdvW9Hr5VajUfT1sGNkHGNsNdc7rnYvxTDrZo/tO3dddu8Sk3GP7aXsO6o1zuOoe7Uel7Ja2z0Ov33n/rznQ/uPm12sJC/7iovyrc+/qOs+lvs5Ney+872OnLtuIv/+uuflqisuXHQ7lqKqlm3G+Y4k26rqsiQfT/LyJN85p8wtSa5P8hdJXpbk/ScKzSyPY2u45r6Ar4ZfZ1ttHKHNqXRsbfVSgsmxNbdZuy7J4g84W41je5j+OmZ6umXy/Gdn6rwLc+cDj3bvP2zda9ZUnn/J01Zs3y5ktbZ7HIY9MHGpY3Oc5nsdOTI5lbsPHF624Nxr6ODcWjtaVa9M8r7MnI7ura21u6rq9Ul2tdZuSXJTkpuram+ST2YmXDNmo3jzoI8jtDnVlhpMhg2+q3VsDxPkjs0aP3blS5OJtXnj++9d1KyxEMmJDDOpNezYHMYwgX2+15H16yZy5ZYNp6q5SzaS8zi31m5Ncuucba+bdfmzSb5tFHUxWl7Al4fZfVaqYYPvmTi2h52lhxMZ5zdISzVsYJ/7OrJ+3US2b92YnZdvPmVtXio/uQ3LwOw+K9WwwXecY3tcX0mvxuUprC7j+gZpqc+pYQP77NeR6dZy5ZYN2Xn55kyswPdIwRmWidl9VqJRBN9xjO1xfiW9WpencPobZmwO85waxYfJY68jL3rW+V3lx2UUP4ACwCp27A3rW59/0bL9wMSwjpvhqjXL+iMmx2bpz167JpUs+49jwEKGGZvDPKeOBfbZTtcPk2acAVh1xrlcwtIrVqphxuYwz6kz6VgHwRmAVWfcyyUsvWKlWurYHOY5dSZ9mBScAVh1zqQZLlgOozhQ+Ez4MCk4A7DqnEkzXLAcPKf6CM4ArEpnygwXLBfPqZNzVg0AAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDJx2pqdbJs9/do5c8tW584FHMz3dxt0kAE4DgjNwyowjwE5Pt7zhvffksStfmiOXfW3e+P5784b33iM8AzA0wRnOAGdSgN2971D2HnwsWbsuqTV54uh09h58LLv3HTql9QJw+hOcYRUYJvieaQH2/kc+k8mj08dtmzw6nfsf+cwprZelsawGWE0EZ1jhhg2+Z1qAvfT8p2Td2uNf2tatXZNLz3/KKa2XxbOsBlhtBGdY4YYNvmdagN2+dWOevfm8nL12TSrJ2WvX5Nmbz8v2rRtPab0snmU1wGqzdtwNON1sWH9WXvSs85e0X5Il7cvp7Y77Pzlv8J1urWu8PD55NH/wVw/m8cmpz29bv24iVz/3i0/peHvBZU/Pn3/0E9m971COTE5l/bqJbN+6MT+081mZWFOnrN4kueVZX5Pb9xzM3QcO58otG7Lz8s2LrtNz8tQbdmwDLDfBGVa452zZkPXrJp4UfK/csqFr/52Xb872rRufFGB3Xr75VDU5STKxpnLzDS8cOsAute6rrrgwV11x4Smvi6UbdmwDLDfBGVa4YYOvAMtKNa4PdQBLJTivAFPTLY9v/JJMPuXC3HbPQ8sWalgdRhF8BVhWonF+qANYCsF5zKamW15x0wfz8LZvSVuzNj/6jg9n+9aNufmGF3rz4PME3+Xlw+zyMbaB1cRZNcbs9j0Hs3vfobSJmaPKH5+cyu59h3L7noPjbhqckWZ/mD100VflR9/x4bzipg9myinSAM54gvOY3XXgcI7MOjAmSY5MTuXuA4fH1CI4s/kwC8BCBOcxO3ZU+WyOKofx8WEWgIUIzmN27Kjyc9dNpJKc66hyGCsfZgFYiIMDx8xR5bCyOEUaAAsRnFcAR5XDyuHDLAALEZwB5vBhFoD5WOMMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADr6r/5HAAAJ2UlEQVQIzgAA0EFwBgCADoIzAAB0GCo4V9XTq+qPqurewd+nzVNme1X9RVXdVVV/WVXfMUydAAAwDsPOOL86yW2ttW1Jbhtcn+vxJN/TWntOkquT/HpVbRyyXgAAWFbDBudrk7x9cPntSV46t0Br7W9ba/cOLh9IcjDJpiHrBQCAZTVscL6wtfZgkgz+bj5R4ap6QZJ1ST66wO03VtWuqtr18MMPD9k0AAAYnbUnK1BV/zPJF89z088upqKqekaSm5Nc31qbnq9Ma+3NSd6cJDt27GiLuX8AADiVThqcW2svXui2qnqoqp7RWntwEIwPLlBuQ5I/SPLa1toHltxaAAAYk2GXatyS5PrB5euT/N7cAlW1LsnvJPmN1tq7h6wPAADGYtjg/MtJvqGq7k3yDYPrqaodVfWWQZlvT/KPknxvVe0e/Ns+ZL0AALCsTrpU40Raa48kuWqe7buSfP/g8m8m+c1h6gEAgHHzy4EAANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAh6GCc1U9var+qKruHfx92gnKbqiqj1fVfximzpVqarrl8Y1fkkPPfFFuu+ehTE23cTcJAIARGnbG+dVJbmutbUty2+D6Qn4xyZ8MWd+KNDXd8oqbPpiHt31LDl30VfnRd3w4r7jpg8IzAMBpZNjgfG2Stw8uvz3JS+crVFVfkeTCJH84ZH0r0u17Dmb3vkNpE+uSWpPHJ6eye9+h3L7n4LibBgDAiAwbnC9srT2YJIO/m+cWqKo1SX4lyU8OWdeKddeBwzkyOXXctiOTU7n7wOExtQgAgFFbe7ICVfU/k3zxPDf9bGcdP5zk1tbavqo6WV03JrkxSS6++OLOux+/52zZkPXrJvL4rPC8ft1ErtyyYYytAgBglE4anFtrL17otqp6qKqe0Vp7sKqekWS+tQkvSvK1VfXDSc5Lsq6qHmutPWk9dGvtzUnenCQ7duxYNQuEd16+Odu3bszufYdyZHIq69dNZPvWjdl5+ZMm4AEAWKVOGpxP4pYk1yf55cHf35tboLX2XccuV9X3JtkxX2hezSbWVG6+4YW5fc/B3H3gcK7csiE7L9+ciTUnnmEHAGD1GDY4/3KSd1XVDUk+luTbkqSqdiT5odba9w95/6vGxJrKVVdcmKuuuHDcTQEA4BSo1lbmiogdO3a0Xbt2jbsZAACc5qrqQ621HScr55cDAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0KFaa+Nuw7yq6uEkD4yp+guSfGJMda9G+mtx9Nfi6K/F0V+Lo78WR38tnj5bnHH11yWttU0nK7Rig/M4VdWu1tqOcbdjtdBfi6O/Fkd/LY7+Whz9tTj6a/H02eKs9P6yVAMAADoIzgAA0EFwnt+bx92AVUZ/LY7+Whz9tTj6a3H01+Lor8XTZ4uzovvLGmcAAOhgxhkAADoIzrNU1dVVtaeq9lbVq8fdntWgqu6vqr+qqt1VtWvc7VlpquqtVXWwqv561ranV9UfVdW9g79PG2cbV5IF+uvnq+rjgzG2u6q+cZxtXEmqamtV/XFV3VNVd1XVjw+2G2PzOEF/GWPzqKpzqur/q6qPDPrrFwbbL6uqDw7G13+rqnXjbutKcIL+eltV/d2s8bV93G1dSapqoqo+XFW/P7i+oseX4DxQVRNJ3pTkmiRXJrmuqq4cb6tWja9rrW1fyaePGaO3Jbl6zrZXJ7mttbYtyW2D68x4W57cX0nya4Mxtr21dusyt2klO5rkVa21K5J8ZZIfGbxuGWPzW6i/EmNsPk8k+frW2pcn2Z7k6qr6yiT/OjP9tS3Jo0luGGMbV5KF+itJfnLW+No9viauSD+e5J5Z11f0+BKcv+AFSfa21u5rrU0meWeSa8fcJla51tqfJvnknM3XJnn74PLbk7x0WRu1gi3QXyygtfZga+3OweVPZ+bN55kxxuZ1gv5iHm3GY4OrZw3+tSRfn+Q9g+3G18AJ+osFVNVFSb4pyVsG1ysrfHwJzl/wzCT7Zl3fHy+oPVqSP6yqD1XVjeNuzCpxYWvtwWTmjTzJ5jG3ZzV4ZVX95WAph2UH86iqS5M8L8kHY4yd1Jz+SoyxeQ2+Rt+d5GCSP0ry0SSHWmtHB0W8V84yt79aa8fG1y8NxtevVdXZY2ziSvPrSX4qyfTg+vlZ4eNLcP6CmmebT4on99WttednZonLj1TVPxp3gzjt/Mckz8rMV58PJvmV8TZn5amq85L89yT/vLV2eNztWenm6S9jbAGttanW2vYkF2Xmm9kr5iu2vK1aueb2V1U9N8nPJPkHSf5hkqcn+ekxNnHFqKpvTnKwtfah2ZvnKbqixpfg/AX7k2yddf2iJAfG1JZVo7V2YPD3YJLfycwLKyf2UFU9I0kGfw+OuT0rWmvtocGb0XSS/xxj7DhVdVZmQuB/ba399mCzMbaA+frLGDu51tqhJLdnZm34xqpaO7jJe+U8ZvXX1YMlQq219kSS/xLj65ivTvKSqro/M8tjvz4zM9ArenwJzl9wR5Jtg6M51yV5eZJbxtymFa2qnlJVTz12Ock/SfLXJ96LzIyr6weXr0/ye2Nsy4p3LAAO/NMYY583WA94U5J7Wmu/OusmY2weC/WXMTa/qtpUVRsHl9cneXFm1oX/cZKXDYoZXwML9NffzPoQW5lZr2t8JWmt/Uxr7aLW2qWZyVzvb619V1b4+PIDKLMMTkH060kmkry1tfZLY27SilZVX5KZWeYkWZvkt/TZ8arqHUl2JrkgyUNJfi7J7yZ5V5KLk3wsybe11hwQlwX7a2dmvkJvSe5P8oPH1u+e6arqa5L8WZK/yhfWCL4mM+t2jbE5TtBf18UYe5Kq+rLMHJw1kZmJtne11l4/eO1/Z2aWHXw4yXcPZlPPaCfor/cn2ZSZZQi7k/zQrIMISVJVO5P8H621b17p40twBgCADpZqAABAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOjw/wOCsl1szxo+SgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 864x576 with 1 Axes>"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAs4AAAHiCAYAAAD4akr8AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzt3X+4XXddJ/r3JydNSCk1hSaFtmmLUmuLj4ZOhh8yOpEfd4qj4PVBBRVxbrVyFX/MMCr+uI7DXB1n5vHqw4XxygwIgwoy6mhH60Ut9ur1QiVARNtObahAYktTKjEtTZPmnO/94+zTnpyeJN+Tvc/e++S8Xs+T5+y19lr7+93fvbL3e33Xd61VrbUAAACntmHSFQAAgLVAcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAOsc1X1yap6yRmu+5VVdeeo6wQwjQRngNOoqluq6nNVtXkF67SqetZq1msSlr6v1tqfttaummSdAMZFcAY4haq6IslXJmlJXj7RypxGVW3smQfAmRGcAU7t25N8KMk7k7x2YeagF/o7F01/R1X9v4PHfzKY/RdV9VBVffNg/ndV1b6q+ruqurGqLl60/rOr6g8Hz91XVT82mL+5qn6hqu4Z/PuFhZ7vqtpdVQeq6keq6jNJfnm5eYNlv7aq9lbVoar6/6rqy5Z7s1X13Kr64GC5e6vqLVW16WTva6G8RetfPWibQ1V1W1W9fNFz76yqt1bV71XVg1V1a1V90Zl9LADjJzgDnNq3J/nVwb9/UlUXnW6F1tpXDR5+eWvtvNbar1fVi5L82yTflOQZST6V5L1JUlVPSfJHSf7vJBcneVaSmwev8eNJnp9kZ5IvT/LcJD+xqLinJ3lqksuT3LDcvKq6Nsk7knx3kqcl+aUkN55k6Mlskn+e5MIkL0jy4iTfc7L3tXjFqjonyX9P8gdJtif5viS/WlWLh3K8Osm/TnJBkn1JfnrZRgSYQoIzwElU1T/KfPh8X2vtI0k+keRbzvDlvjXJO1prH22tHU3yo0leMBgK8rVJPtNa+7nW2iOttQdba7cuWu9NrbWDrbX7Mx86X7PodeeS/KvW2tHW2pGTzPuuJL/UWru1tTbbWntXkqOZD+QnaK19pLX2odba8dbaJzMfsv9x53t8fpLzkvxsa+1Ya+0DSX4382F5wW+11v68tXY88zsjOztfG2DiBGeAk3ttkj9orX12MP1rWTRcY4Uuznwvc5KktfZQkgeSXJJkR+ZD+WnXGzy+eNH0/a21R5ass3Te5UneMBg+caiqDg3KvHjJeqmqL66q362qz1TV4SQ/k/ne5x4XJ9nfWptbUt9LFk1/ZtHjhzMftAHWBCeNACyjqrZkfljFzGCscJJsTrK1qr48yeeTnLtolaef5iXvyXyAXXj9J2d+2MTfJtmfE3tll1vvtsH0ZYN5C9oy6yydtz/JT7fWeoZF/GKSjyV5dWvtwar6wSSv7Fhvoa47qmrDovB8WZK/7lwfYKrpcQZY3tdnfrzvNZkfTrAzydVJ/jTz4573JvmGqjp3cHm265esf1+SL1w0/WtJ/llV7RyMLf6ZJLcOhkP8bpKnV9UPDk4GfEpVPW+w3nuS/ERVbauqC5P8ZJJfWeF7+U9JXldVz6t5T66qfzoYW73UU5IcTvJQVX1Jkv/1NO9rsVszv0Pxw1V1TlXtTvJ1GYzlBljrBGeA5b02yS+31j7dWvvMwr8kb8n8uOOfT3Is80HyXZkfr7vYTyV512BoxDe11m5O8r8l+c0k9yb5oiSvSpLW2oNJXpr5kPmZJHcl+erB6/zvSfYk+XiSv0zy0cG8bq21PZkf5/yWJJ/L/El533GSxf9l5sdxP5j5wP3rS54/4X0tKedY5i/Z97Ikn03yH5N8e2vtf6ykvgDTqlpb7igfAACwmB5nAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADlN7A5QLL7ywXXHFFZOuBgAAZ7mPfOQjn22tbTvdclMbnK+44ors2bNn0tUAAOAsV1Wf6lnOUA0AAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0GElwrqp3VNXBqvqrkzxfVfXmqtpXVR+vqmtHUS4AAIzLqHqc35nkulM8/7IkVw7+3ZDkF0dU7kjNzrXcfMd9efPNd+XmO+7L7FybdJUAAJgSG0fxIq21P6mqK06xyCuS/JfWWkvyoaraWlXPaK3dO4ryR2F2ruU1b781e/cfypFjs9myaSY7d2zNu69/XmY21KSrBwDAhI1rjPMlSfYvmj4wmDc1brnzYPbuP5SHj82mJXn42Gz27j+UW+48OOmqAQAwBcYVnJfrsn3COIiquqGq9lTVnvvvv38M1XrcbfcczpFjsyfMO3JsNrffc3is9QAAYDqNKzgfSLJj0fSlSe5ZulBr7W2ttV2ttV3btm0bU9XmPfvi87Nl08wJ87Zsmsk1F58/1noAADCdxhWcb0zy7YOrazw/yd9P0/jmJNl91fbs3LE1NXssaXM5dzDGefdV2yddNQAApsBITg6sqvck2Z3kwqo6kORfJTknSVpr/1eSm5J8TZJ9SR5O8s9GUe4ozWyovPv65+UF33B9jj15e37uJ/55dl+13YmBAAAkGd1VNV59mudbku8dRVmraWZD5dxDd+fcQ3fnxVdfNOnqAAAwRdw5EAAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQIeRBOequq6q7qyqfVX1xmWev6yq/riqPlZVH6+qrxlFuQAAMC5DB+eqmkny1iQvS3JNkldX1TVLFvuJJO9rrT0nyauS/MdhywUAgHEaRY/zc5Psa63d3Vo7luS9SV6xZJmW5PzB4y9Ics8IygUAgLEZRXC+JMn+RdMHBvMW+6kk31ZVB5LclOT7lnuhqrqhqvZU1Z77779/BFUDAIDRGEVwrmXmtSXTr07yztbapUm+Jsm7q+oJZbfW3tZa29Va27Vt27YRVA0AAEZjFMH5QJIdi6YvzROHYlyf5H1J0lr7YJInJblwBGUDAMBYjCI4fzjJlVX1zKralPmT/25cssynk7w4Sarq6swHZ2MxAABYM4YOzq2140len+T9Se7I/NUzbquqN1XVyweLvSHJd1XVXyR5T5LvaK0tHc4BAABTa+MoXqS1dlPmT/pbPO8nFz2+PckLR1EWAABMgjsHAgBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB1GEpyr6rqqurOq9lXVG0+yzDdV1e1VdVtV/dooygUAgHHZOOwLVNVMkrcmeWmSA0k+XFU3ttZuX7TMlUl+NMkLW2ufq6rtw5YLAADjNIoe5+cm2ddau7u1dizJe5O8Ysky35Xkra21zyVJa+3gCMoFAICxGUVwviTJ/kXTBwbzFvviJF9cVX9WVR+qqutGUC4AAIzN0EM1ktQy89oy5VyZZHeSS5P8aVV9aWvt0AkvVHVDkhuS5LLLLhtB1QAAYDRG0eN8IMmORdOXJrlnmWV+p7X2aGvtb5LcmfkgfYLW2ttaa7taa7u2bds2gqoBAMBojCI4fzjJlVX1zKralORVSW5cssxvJ/nqJKmqCzM/dOPuEZQNAABjMXRwbq0dT/L6JO9PckeS97XWbquqN1XVyweLvT/JA1V1e5I/TvJDrbUHhi0bAADGZRRjnNNauynJTUvm/eSixy3Jvxj8AwCANcedAwEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQYSXCuquuq6s6q2ldVbzzFcq+sqlZVu0ZRLgAAjMvQwbmqZpK8NcnLklyT5NVVdc0yyz0lyfcnuXXYMgEAYNxG0eP83CT7Wmt3t9aOJXlvklcss9y/SfLvkzwygjIBAGCsRhGcL0myf9H0gcG8x1TVc5LsaK397gjKAwCAsRtFcK5l5rXHnqzakOTnk7zhtC9UdUNV7amqPffff/8IqgYAAKMxiuB8IMmORdOXJrln0fRTknxpkluq6pNJnp/kxuVOEGytva21tqu1tmvbtm0jqBoAAIzGKILzh5NcWVXPrKpNSV6V5MaFJ1trf99au7C1dkVr7YokH0ry8tbanhGUDQAAYzF0cG6tHU/y+iTvT3JHkve11m6rqjdV1cuHfX0AAJgGG0fxIq21m5LctGTeT55k2d2jKBMAAMbJnQMBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0GElwrqrrqurOqtpXVW9c5vl/UVW3V9XHq+rmqrp8FOUCAMC4DB2cq2omyVuTvCzJNUleXVXXLFnsY0l2tda+LMlvJPn3w5YLAADjNIoe5+cm2ddau7u1dizJe5O8YvECrbU/bq09PJj8UJJLR1AuAACMzSiC8yVJ9i+aPjCYdzLXJ/n95Z6oqhuqak9V7bn//vtHUDUAABiNUQTnWmZeW3bBqm9LsivJf1ju+dba21pru1pru7Zt2zaCqgEAwGhsHMFrHEiyY9H0pUnuWbpQVb0kyY8n+cettaMjKBcAAMZmFD3OH05yZVU9s6o2JXlVkhsXL1BVz0nyS0le3lo7OIIyAQBgrIYOzq2140len+T9Se5I8r7W2m1V9aaqevlgsf+Q5Lwk/7Wq9lbVjSd5OQAAmEqjGKqR1tpNSW5aMu8nFz1+ySjKAQCASXHnQAAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdRnIDFNan2bmWW+48mNvuOZxnX3x+dl+1PTMbatLVAgBYFYIzZ2R2ruU1b781e/cfypFjs9myaSY7d2zNu69/nvAMAJyVDNXgjNxy58Hs3X8oDx+bTUvy8LHZ7N1/KLfceXDSVQMAWBWCM2fktnsO58ix2RPmHTk2m9vvOTyhGgEArC7BmTPy7IvPz5ZNMyfM27JpJtdcfP6EagQAsLoEZ87I7qu2Z+eOranZY0mby7mDMc67r9o+6aoBAKwKwZkzMrOh8u7rn5dtd/33bD3wZ/k/X/0cJwYCAGc1V9XgjM1sqJx76O6ce+juvPjqiyZdHQCAVSU4A1PJdcIBmDaCMzB11ut1wltraS1pi6Yfe+6xeYuWH8xdPG+chi23ZfkXWPq6S5dqSxY4VTW66jim9jvZ+13rRrn9jauFlm5Dy5V9qvd1ys/yNG+ia5McQ0NM4/Z47jkb8wXnnjPpapyS4AxMhYXQONdabr7j8euEJ49fJ/ymv7w3X/XF25I2/6W/EDJba5kbzMtj8+an59qiQDp4/ZbB35bksWUenze33I9q52/M8bm5fPRTn8snPvv5fOGFT85zLrsgG+rxOuWEeg/qtoLXBzhbbT9/s+DM+jA71x4LL4vDx1x7PLgs+zftNHv1q2e5HofVK2tsRa2a1p4YVpNTBMI8HlYXHs8t2UaWbjML/uiO+5a9Tvif/PX9ufC8zWN6xys3N9fyM79/R/YdfCjHjs9l08YNedb28/JjL7s6G87innKA9UJwnrDjs3N58JHj+fsjj+bzx46feBh28Hjx4ZQTD9MuzHvi4dyly+Y0r7H8Oqc/JHr4yKNJkj//m78LZ6e5uZa9+w/lkw98Plc87cnZuWPrqofAK5725GzauCFHj889Nm/Txg254mlPXtVyh7V3/6HsO/jQY/U+enwu+w4+lL37D+Xayy+YcO0AGJbgPGazcy0PPvJoDh85nsOPPJqHjh4/K3ojOTtNqgd1546tedb283Lbpz+bzGzM5nM25lnbz8vOHVtXrcxR+OQDn8+xRWE/SY4dn8snH/i84AxwFhCcV9ncXMuDR4/n8JFH53uVjx7PnKDMGjGpHtQNGyo/9rKr890/8IbMnndRXv+6G8bS0z2stdpTDkAfwXnEWmt56Oj80IvDR47nwUceFZRZsybZg7phQ2XTA/uSB/bl2st/ZFXLGpW12lMOQB/BecT27j+URx6dO/2CsAboQV2ZtdpTDkAft9weMb3LnE0WelBz/FjS5rJ5MMZZD+rJLfSUb/nUn+Xayy8QmgHOIoIzcFILPajn3f7b2fI3f5rvf9GVLq0GwLolOAOnpAcVAOYJzgAA0EFwBgCADoIzAAB0cDm6KTDsLY0ncUtkAID1RnCesGFvaTypWyIDAKw3hmpM2OJbGreceEvjcawPAEAfwXnCTnVL43GsDwBnq7m5lo9+6nP5rY8eyEc/9bnMuUsZQzJUY8KGvaWxWyIDwBMZyshq0OM8YcPe0tgtkQHgiQxlZDUIzhM27C2N3RIZAJ5o0kMZDRM5OxmqMQUWbmmcB/bl2st/ZOzrT4JL6AHD8j3CqUxyKKNhImcvwZmx84UCDMv3CKezMJTxtk9/NpnZmM3nbBzbUMbFw0SSE4eJXHv5BatePqvHUA3GzrgzYFi+R8ZrLQ47mORQxkkPE2H16HFm7E71hTLte+IODcN0WMvfI2vNWu7dn9RQRle8OnsJzozdWv1CWcs/HnC2WavfI2uRYQcrN8lhIjp4VpfgzNhN8gslOfMvFT8eMD0m/T2ynujdX7mFYSLf/QNvyOx5F+X1r7thLAFWB8/qE5wZu0l9oSTDfan48YDpMcnvkfVG7/6ZmcQwkUl28KyXnm7BmYmY1LizYb5U/HjAdFmLl+Jci/Turx2T6uAZRU/33FzLBz/xQH7v4/fm2Refn91Xbc/MFAZvwZl1ZZgvlbX847FeegKA0dO7v3ZMqoNn2J7uheB99/0P5ZFH57Jl00x27tiad1//vKkLzyO5HF1VXVdVd1bVvqp64zLPb66qXx88f2tVXTGKcmGlFr5UFuv9Ulmrd2lc+EJ68wfuym985EDe/IG78jO/f8eauJwUMB0Weve3fOrPcu3lF4z1e28tXgpvUhY6eHL8WNLmsnnQ89vbwXOmbT3s5fcWgveRR+cvL/nwsdns3X8ot9x5sGv9cRq6x7mqZpK8NclLkxxI8uGqurG1dvuixa5P8rnW2rOq6lVJ/l2Sbx62bFipYXuNhzk0PKleXyc1rg+OKqwdPqt+TnZbmWGODgzT1sP2dC8XvI8cm83t9xzOi6++qOs1xmUUQzWem2Rfa+3uJKmq9yZ5RZLFwfkVSX5q8Pg3krylqqq1ZreRFRvmR2c9nunspMazn3CxdvisVsaO/8qdaQfPMG09bKfUcsF7y6aZXHPx+d31H5caNrtW1SuTXNda+87B9GuSPK+19vpFy/zVYJkDg+lPDJb57Mle96mXX91e+mPvGKpuZ2LvX+xNkuz88p1ntP5DR49n7gza9K7b/ypJcuU1X3pG5Q67/pkaptwzWbe1lk//3ZEceXQ2rSVVyZZzZnLZU7ekqv9HZ9z1fvCR4/nbQ0eyeNOoSi7ZuiVPedLqnmowirLX4vY1SeOu9yS3r0lba9vIWv6sJtHW9z94NJ996NgT5m87b1MufMrmrteY5Daylr47h23r1lr+et/dycymXHzxM3Le5pnu3+WF3/aHjz6apLJhQ+W8zRvzJU9/yop+24fxvtd9xUdaa7tOt9wo/pcu946WJseeZVJVNyS5IUnOe8YXDV+zM3CmgTmZD92zcy1XXr3y/yDD/qcaZv1h/mMPU+6ZrPvQ0dnHQnOStJYceXQ2Dx2dXdGPzrjr/ciiOi9oLTn6aH+9z/RzOm/zTLacM/OEnY3zNs90v8Za3L4muTM67npPcvsaxfprqa2HXX/Yz2qYcHKmdV4w7O/UmZT9pHNmUpUn7GhsPmf6v78mWfaZrDNsW1dVrrryzLJbVeWyp27J0Uc359HZlnM3z2TrlnPGFppXYhQ9zi9I8lOttX8ymP7RJGmt/dtFy7x/sMwHq2pjks8k2XaqoRq7du1qe/bsGapu47Z79+48+MjxvOVXf2fSVVmR7/2WlydJ3vprN064Jqf3Wx89kN/4yIET9roqySv/waX5hmsvnVS1Tuujn/pc3vyBu044DLV544Z8/4uu7D7cOMznNMkxlZPavoYtdy3Ve9Lb17DrT6Kt5+ba0EO2xv1ZLQzzWHo4vHeYxyje8zDOpL0mPbRlkr+P4y570m2dJNvP35wv2nbeWMpaqqrG1uP84SRXVtUzk/xtklcl+ZYly9yY5LVJPpjklUk+YHwzZ2KtXkt5YfzX0i+kcV3KbsOGyrWXX2BM4Flq0tvXWrMQEB665uuTmY158wfuGltAGOazWhiDmo2bkqxsDOok3/NC+cee9qzMnndRPvqpz3WH9oXzUtbbyZRn2l7DWK9tvVJDB+fW2vGqen2S9yeZSfKO1tptVfWmJHtaazcmeXuSd1fVviR/l/lwDSu2VgOCLyRWk+1rZYYJoMMa5rMa5kTfSb7nYUP7etvxn+ROznpr6zMxkjMRWms3JblpybyfXPT4kSTfOIqyWN/WckDwhcRqGmb7mkTv1iRN+kozZ/pZDXPEbZLveZKhfS3SXtNtJDdAgXFa+NH5hmsvHfuF+OFss7h368gzv/KMbpCzELyPXP7CNXGDimFuhLRgEu954Yjb5o0bUsmKbm4xivd8poa9OcZ6o72m23Rf+wZIsv56BCdtPbX3sL1bkx47eyaGHfI1qfc8zBG3SQ5zW6vnpkyK9ppugjNMubUYTNay9dbewx7CX4uHlYcd8jXpMdJnMsxjksPc1uq5KZPagV6r7bVeCM4w5dZiMFnL1lt7r8atctfCnSmHGRO+Ht/zsOWutXNTJn2C3lprr/XEGGeYcsa7jdd6a+9hxs0mkx07Oynr8T0Pa62dm3LCDnRtOGEHehzWWnutJ3qcYcoZ7zZe6629h+3dWo+Hldfje15v1upRBVaf4LzOraeToNYqP9LjtR7be5hD+OvxsPJ6fM/rzXrbgaaf4LyOrbeToNYqP9Ljpb1Xbj1eQ9p12c9u63EHmj6C8zq23k6CWsv8SI+X9h4PO+9MKzvQnIyTA9ex9XYSFDBdJn0CFpyKE/RYjuC8jjkzHJgkO+/AWiM4r2PDXoYKGK21duvqYdl5B9YaY5zXMWO4YHqsx/G+TsAC1hrBeZ1zEhRMh/V4sq6dd2CtEZwBpsB6veGCnXdgLTHGGWAKGO8LMP0EZ4Ap4GRdgOlnqAbAFDDeF2D6Cc4AU8J4X4DpZqgGAAB0EJwBAKCD4DxiGzdUZoxJBAA46xjjPGJbNs3kH15xQR48ejyHjzyaw0eO58FHHs1ZfudcAICznuC8Cqoq5z/pnJz/pHOSC+ZvpbsQpP/+yKN56OjxNEEaAGBNEZzHYMOGyhdsOSdfsOWc7EgyO9fy4CPzvdF/f+TRfP6YIA0AMO0E5wmY2VDZeu6mbD13U5Lk+OxcDj9yPJ9f0hPdMj9x4rxFjwdPnDhvcUlt2fknX/7xMk/6/DLLtzY/v6VlriVzraW1weO5ZpgKAHBWEJynwMaZDXnqkzflqU/eNOmqrJr5AD0fy+eD9eN/T5Wr22m64tdKJh/7EYVVKO+xnaTMfy4LRSzsNOWx505cdvM5M0mSS7ZuyVxrg3/zr7iwozU3N/83yWPPP7Z9DMpavM0AwCQIzozFhg2V+RsJs95s3jh/8Z7LnnbuyF6ztSeG+MWh+vhsy8zTr8rRJ2/PAw8dzQuvvDAbUo8dFVk4IrIQ+BceLzy3cBTlhDI790aWHt1ZHPxzws7HE4/YLKz72N9ljhqt7KjTahuusJPVdensJx75WnJk7AzKADgTgjOw5lRV6rH9sBN3yGbnWv6Xd96az37x16Vt2Jgf/s2PZ+eOrXn39c9zqUiWdbojW6MrZyzFjN0o39bYPotly166zMnrMkw1e9bt3VEfxjRujxtnpv87WnAGziq33Hkwe/cfSpuZH/r08LHZ7N1/KLfceTAvvvqiCdeOaVQ1nh/rMRWzxmkkppsboABnldvuOZwjx2ZPmHfk2Gxuv+fwhGoEwNlVV5QqAAAKy0lEQVRCcAZWzexcy8NbvzCHLnlBbr7jvsyO4RIrz774/GzZNHPCvC2bZnLNxeevetkAnN0EZ2BVzM61vObtt+b+K78uhy79inzfez6W17z91lUPz7uv2p6dO7bm3E0zqSTnbprJzh1bs/uq7ataLgBnP2OcgVUxqbHGMxsq777+ebnlzoO5/Z7Duebi87P7qu1ODARgaIIzsCpONdZ4tU/Sm9lQefHVFzkZEICRMlQDWBXGGgNwthGcgVVhrDEAZxtDNYBVYawxAGcbwRlYNcYaA3A2MVQDAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgwVnKvqqVX1h1V11+DvBcsss7OqPlhVt1XVx6vqm4cpEwAAJmHYHuc3Jrm5tXZlkpsH00s9nOTbW2vPTnJdkl+oqq1DlgsAAGM1bHB+RZJ3DR6/K8nXL12gtfbXrbW7Bo/vSXIwybYhywUAgLEaNjhf1Fq7N0kGf7efauGqem6STUk+MWS5AAAwVhtPt0BV/VGSpy/z1I+vpKCqekaSdyd5bWtt7iTL3JDkhiS57LLLVvLyAACwqk4bnFtrLznZc1V1X1U9o7V27yAYHzzJcucn+b0kP9Fa+9Apynpbkrclya5du9rp6gYAAOMy7FCNG5O8dvD4tUl+Z+kCVbUpyX9L8l9aa/91yPIAAGAihg3OP5vkpVV1V5KXDqZTVbuq6j8PlvmmJF+V5Duqau/g384hywUAgLE67VCNU2mtPZDkxcvM35PkOwePfyXJrwxTDgAATJo7BwIAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdhgrOVfXUqvrDqrpr8PeCUyx7flX9bVW9ZZgyp9XsXMvDW78why55QW6+477MzrVJVwkAgBEatsf5jUlubq1dmeTmwfTJ/Jsk/8+Q5U2l2bmW17z91tx/5dfl0KVfke97z8fymrffKjwDAJxFhg3Or0jyrsHjdyX5+uUWqqp/kOSiJH8wZHlT6ZY7D2bv/kNpM5uS2pCHj81m7/5DueXOg5OuGgAAIzJscL6otXZvkgz+bl+6QFVtSPJzSX5oyLKm1m33HM6RY7MnzDtybDa333N4QjUCAGDUNp5ugar6oyRPX+apH+8s43uS3NRa219VpyvrhiQ3JMlll13W+fKT9+yLz8+WTTN5eFF43rJpJtdcfP4EawUAwCidNji31l5ysueq6r6qekZr7d6qekaS5cYmvCDJV1bV9yQ5L8mmqnqotfaE8dCttbcleVuS7Nq1a80MEN591fbs3LE1e/cfypFjs9myaSY7d2zN7que0AEPAMAaddrgfBo3Jnltkp8d/P2dpQu01r514XFVfUeSXcuF5rVsZkPl3dc/L7fceTC333M411x8fnZftT0zG07dww4AwNoxbHD+2STvq6rrk3w6yTcmSVXtSvK61tp3Dvn6a8bMhsqLr74oL776oklXBQCAVVCtTeeIiF27drU9e/ZMuhoAAJzlquojrbVdp1vOnQMBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQoVprk67Dsqrq/iSfmlDxFyb57ITKXou018por5XRXiujvVZGe62M9lo5bbYyk2qvy1tr20630NQG50mqqj2ttV2Trsdaob1WRnutjPZaGe21MtprZbTXymmzlZn29jJUAwAAOgjOAADQQXBe3tsmXYE1RnutjPZaGe21MtprZbTXymivldNmKzPV7WWMMwAAdNDjDAAAHQTnRarquqq6s6r2VdUbJ12ftaCqPllVf1lVe6tqz6TrM22q6h1VdbCq/mrRvKdW1R9W1V2DvxdMso7T5CTt9VNV9beDbWxvVX3NJOs4TapqR1X9cVXdUVW3VdUPDObbxpZxivayjS2jqp5UVX9eVX8xaK9/PZj/zKq6dbB9/XpVbZp0XafBKdrrnVX1N4u2r52Trus0qaqZqvpYVf3uYHqqty/BeaCqZpK8NcnLklyT5NVVdc1ka7VmfHVrbec0Xz5mgt6Z5Lol896Y5ObW2pVJbh5MM++deWJ7JcnPD7axna21m8Zcp2l2PMkbWmtXJ3l+ku8dfG/ZxpZ3svZKbGPLOZrkRa21L0+yM8l1VfX8JP8u8+11ZZLPJbl+gnWcJidrryT5oUXb197JVXEq/UCSOxZNT/X2JTg/7rlJ9rXW7m6tHUvy3iSvmHCdWONaa3+S5O+WzH5FkncNHr8rydePtVJT7CTtxUm01u5trX108PjBzP/4XBLb2LJO0V4so817aDB5zuBfS/KiJL8xmG/7GjhFe3ESVXVpkn+a5D8PpitTvn0Jzo+7JMn+RdMH4gu1R0vyB1X1kaq6YdKVWSMuaq3dm8z/kCfZPuH6rAWvr6qPD4ZyGHawjKq6Islzktwa29hpLWmvxDa2rMFh9L1JDib5wySfSHKotXZ8sIjfykWWtldrbWH7+unB9vXzVbV5glWcNr+Q5IeTzA2mn5Yp374E58fVMvPsKZ7eC1tr12Z+iMv3VtVXTbpCnHV+MckXZf7Q571Jfm6y1Zk+VXVekt9M8oOttcOTrs+0W6a9bGMn0Vqbba3tTHJp5o/MXr3cYuOt1fRa2l5V9aVJfjTJlyT5h0memuRHJljFqVFVX5vkYGvtI4tnL7PoVG1fgvPjDiTZsWj60iT3TKgua0Zr7Z7B34NJ/lvmv1g5tfuq6hlJMvh7cML1mWqttfsGP0ZzSf5TbGMnqKpzMh8Cf7W19luD2baxk1iuvWxjp9daO5TklsyPDd9aVRsHT/mtXMai9rpuMESotdaOJvnl2L4WvDDJy6vqk5kfHvuizPdAT/X2JTg/7sNJrhyczbkpyauS3DjhOk21qnpyVT1l4XGS/ynJX516LTK/Xb128Pi1SX5ngnWZegsBcOB/jm3sMYPxgG9Pckdr7f9Y9JRtbBknay/b2PKqaltVbR083pLkJZkfF/7HSV45WMz2NXCS9vofi3ZiK/PjdW1fSVprP9pau7S1dkXmM9cHWmvfminfvtwAZZHBJYh+IclMkne01n56wlWaalX1hZnvZU6SjUl+TZudqKrek2R3kguT3JfkXyX57STvS3JZkk8n+cbWmhPictL22p35Q+gtySeTfPfC+N31rqr+UZI/TfKXeXyM4I9lftyubWyJU7TXq2Mbe4Kq+rLMn5w1k/mOtve11t40+O5/b+aHHXwsybcNelPXtVO01weSbMv8MIS9SV636CRCklTV7iT/srX2tdO+fQnOAADQwVANAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB3+f6iy/qYb4GkIAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 864x576 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAs4AAAHiCAYAAAD4akr8AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzt3X+43VddJ/r3JydNm1JioE0qoekPIXZaeDRgBsQfM9Hi3NYflHFQqYrVW60+ijojV0XkQcWRcWauP4YZ7txhKAPWEQYYf3S0XHSKVccr3IYS1LbGhtqSkNqU0hBKQw85Z90/zg6cnJ6TrHP2ztnnJK/X8+Q5e3/3+u619srae7/32uv73dVaCwAAcGJrxt0AAABYDQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzQIeqek1VvaWz7Nuq6l+e6jatdFX1vVX1v4bY/71Vdf0o2wQwDMEZOC1U1f1VdaSqHquqh6rqv1TVeUu8r51VtX/2ttbaG1pr3z+a1n6+jlZVP7XI/X6+qn5zVO1YKeZ7XK21a1prbx9XmwDmEpyB08m3tNbOS/L8JP8wyWsXewdVtXbkrZrf9Uk+Ofi7otWMNSfbBnC686IHnHZaax9P8t4kz02Sqvq+qrqnqj5dVfdV1Q8eK3tsdrmqfrqq/j7JOwb7bhnMXj9WVVvmzohW1bur6u+r6lNV9adV9Zze9lXVuUleluRHkmyrqh1z2zOn/P1V9eKqujrJa5J8x6BdHxncvqWqbqmqT1bV3qr6gVn7TgyWmXx08Pg/VFVbB7d9VVXdMXgMd1TVV83a7/aq+qWq+vMkjyf5kgW2fVFV3VRVD1bVx6vqX1bVxAKP+99V1b6qOjxox9cOti/0uG6vqu8fXF5TVa+tqgeq6mBV/UZVfdHgtksHs/fXV9XHquoTVfWzvf8fAL0EZ+C0MwiG35jkw4NNB5N8c5INSb4vya9V1fNn7fLFSZ6e5JIk35PkmiQHWmvnDf4dmKea9ybZlmRzkjuT/NdFNPGfJXksybuTvG9Q50m11v6fJG9I8t8G7frywU3vSLI/yZbMBPI3VNVVg9t+Isl1memPDUn+9ySPV9XTk/xBkjcmOT/Jryb5g6o6f1aVr0hyY5KnJnlggW1vT3I0ybOTPC/JP0my0JKWO5Jsz0xf/1aSd1fVOSd4XLN97+Df1yX5kiTnJfkPc8p8TZLLk1yV5HVVdcUC7QBYEsEZOJ38blUdSvK/kvxJZsJYWmt/0Fr7aJvxJ0n+MMnXztpvOsnPtdaeaK0d6amotfbW1tqnW2tPJPn5JF9+bAa0w/WZCYlTmQmQ11XVWZ37HmfwIeFrkvx0a+2zrbXdSd6SmYCbzITY17bW9gwe/0daa48k+aYk97bWbm6tHW2tvSPJ3yT5lll3/7bW2l2D2z83d1tmAvA1Sf55a+0zrbWDSX4tycvna2tr7Tdba48M7u9XkpydmaDb47uS/Gpr7b7W2mNJfibJy+csrfmF1tqR1tpHknwkyXwBHGDJBGfgdPLS1trG1tolrbUfPhaCq+qaqvrAYCnDoczMvl4wa7+HW2uf7a1ksPzhlwfLHw4nuX9w0wUn2O3YvlszM2t6bIb695Kck5kguxRbknyytfbpWdseSPLMweWtST66wH4PzNk2e78k2TfPfrO3XZLkrCQPVtWhQd/+p8zMwj9JVb1qsGTmU4OyX5SOPlugvQ8kWZvkwlnb/n7W5cczMysNMDKCM3Baq6qzk/z3JP9nkgtbaxuT3JqkZhVrc3abe32u70xybZIXZyb8XXqsuo4mvSIzr73/Y7Cm+r7MBOdjyzU+k+TcWe2fSLLpBG07kOTpVfXUWdsuTvLxweV9SZ41TzsOZCb4zjZ7v/nqmrttX5Inklww+MCysbW2obX2pPXeg/XMP53k25M8bfD/8Kl8oc9O1udz23txZpaIPHSS/QBGRnAGTnfrMrMk4OEkR6vqmsyswz2Rh5Kcf4KlF0/NTGB8JDMh9w2LaM/3JPmFzKz1PfbvnyX5psH64r9Nck5VfdNg+cZrB+2f3bZLj53RorW2L8n/m+RfVdU5VfVlSW7IF2a035LkF6tq2+BMGF82qOfWJF9aVd9ZVWur6juSXJnk93sfSGvtwcwse/mVqtowOIDvWVX1j+cp/tTMBN2Hk6ytqtdlZs31vI9rHu9I8i+q6rKaOc3gsTXRR3vbCzAswRk4rQ2WMPxYkncleTQzs8W3nGSfv8lMULtvsARhy5wiv5GZpQIfT3J3kg/0tKWqvjIzs9Nvaq39/ax/tyTZm+S61tqnkvxwZgLvxzMzAz37LBvvHvx9pKruHFy+bnC/B5L8TmbWa//R4LZfHTz2P0xyOMlNSdYP1jl/c5JXZeYDwE8l+ebW2id6Hsss35OZDyd3Z6Z/35PkGfOUe19mDqj828z03Wdz/LKP+R7XbG9NcnOSP03yd4P9f3SRbQUYSrV2sm/HAAAAM84AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAd1p68yHhccMEF7dJLLx13MwAAOM196EMf+kRrbdPJyq3Y4HzppZdm165d424GAACnuap6oKecpRoAANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoMJLgXFVvraqDVfXXC9xeVfXGqtpbVX9ZVc8fRb0AALBcRjXj/LYkV5/g9muSbBv8uzHJfxxRvSM1Nd1y2z0P5Y233Zvb7nkoU9Nt3E0CAGCFWDuKO2mt/WlVXXqCItcm+Y3WWkvygaraWFXPaK09OIr6R2FquuUVN30wu/cdypHJqaxfN5HtWzfm5htemIk1Ne7mAQAwZsu1xvmZSfbNur5/sG3FuH3PwezedyiPT06lJXl8ciq79x3K7XsOjrtpAACsAMsVnOebsn3SOoiqurGqdlXVrocffngZmvUFdx04nCOTU8dtOzI5lbsPHF7WdgAAsDItV3Den2TrrOsXJTkwt1Br7c2ttR2ttR2bNm1apqbNeM6WDVm/buK4bevXTeTKLRuWtR0AAKxMyxWcb0nyPYOza3xlkk+tpPXNSbLz8s3ZvnVjamoyadM5d7DGeeflm8fdNAAAVoCRHBxYVe9IsjPJBVW1P8nPJTkrSVpr/3eSW5N8Y5K9SR5P8n2jqHeUJtZUbr7hhXnRt96Qyadszq+89l9k5+WbHRgIAECS0Z1V47qT3N6S/Mgo6jqVJtZUzj10X849dF+uuuLCcTcHAIAVxC8HAgBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6DCS4FxVV1fVnqraW1Wvnuf2i6vqj6vqw1X1l1X1jaOoFwAAlsvQwbmqJpK8Kck1Sa5Mcl1VXTmn2GuTvKu19rwkL0/yfw1bLwAALKdRzDi/IMne1tp9rbXJJO9Mcu2cMi3JhsHlL0pyYAT1AgDAshlFcH5mkn2zru8fbJvt55N8d1XtT3Jrkh+d746q6saq2lVVux5++OERNA0AAEZjFMG55tnW5ly/LsnbWmsXJfnGJDdX1ZPqbq29ubW2o7W2Y9OmTSNoGgAAjMYogvP+JFtnXb8oT16KcUOSdyVJa+0vkpyT5IIR1A0AAMtiFMH5jiTbquqyqlqXmYP/bplT5mNJrkqSqroiM8HZWgwAAFaNoYNza+1oklcmeV+SezJz9oy7qur1VfWSQbFXJfmBqvpIknck+d7W2tzlHAAAsGKtHcWdtNZuzcxBf7O3vW7W5buTfPUo6gIAgHHwy4EAANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAh5EE56q6uqr2VNXeqnr1AmW+varurqq7quq3RlEvAAAsl7XD3kFVTSR5U5JvSLI/yR1VdUtr7e5ZZbYl+ZkkX91ae7SqNg9bLwAALKdRzDi/IMne1tp9rbXJJO9Mcu2cMj+Q5E2ttUeTpLV2cAT1AgDAshlFcH5mkn2zru8fbJvtS5N8aVX9eVV9oKquHkG9AACwbIZeqpGk5tnW5qlnW5KdSS5K8mdV9dzW2qHj7qjqxiQ3JsnFF188gqYBAMBojGLGeX+SrbOuX5TkwDxlfq+19rnW2t8l2ZOZIH2c1tqbW2s7Wms7Nm3aNIKmAQDAaIwiON+RZFtVXVZV65K8PMktc8r8bpKvS5KquiAzSzfuG0HdAACwLIYOzq21o0lemeR9Se5J8q7W2l1V9fqqesmg2PuSPFJVdyf54yQ/2Vp7ZNi6AQBguYxijXNaa7cmuXXOttfNutyS/MTgHwAArDp+ORAAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAECHkQTnqrq6qvZU1d6qevUJyr2sqlpV7RhFvQAAsFyGDs5VNZHkTUmuSXJlkuuq6sp5yj01yY8l+eCwdQIAwHIbxYzzC5Lsba3d11qbTPLOJNfOU+4Xk/ybJJ8dQZ0AALCsRhGcn5lk36zr+wfbPq+qnpdka2vt9090R1V1Y1XtqqpdDz/88AiaBgAAozGK4FzzbGufv7FqTZJfS/Kqk91Ra+3NrbUdrbUdmzZtGkHTAABgNEYRnPcn2Trr+kVJDsy6/tQkz01ye1Xdn+Qrk9ziAEEAAFaTUQTnO5Jsq6rLqmpdkpcnueXYja21T7XWLmitXdpauzTJB5K8pLW2awR1AwDAshg6OLfWjiZ5ZZL3Jbknybtaa3dV1eur6iXD3j8AAKwEa0dxJ621W5PcOmfb6xYou3MUdQIAwHLyy4EAANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAh5EE56q6uqr2VNXeqnr1PLf/RFXdXVV/WVW3VdUlo6gXAACWy9DBuaomkrwpyTVJrkxyXVVdOafYh5PsaK19WZL3JPk3w9YLAADLaRQzzi9Isre1dl9rbTLJO5NcO7tAa+2PW2uPD65+IMlFI6gXAACWzSiC8zOT7Jt1ff9g20JuSPLeEdQLAADLZu0I7qPm2dbmLVj13Ul2JPnHC9x+Y5Ibk+Tiiy8eQdMAAGA0RjHjvD/J1lnXL0pyYG6hqnpxkp9N8pLW2hPz3VFr7c2ttR2ttR2bNm0aQdMAAGA0RhGc70iyraouq6p1SV6e5JbZBarqeUn+U2ZC88ER1AkAAMtq6ODcWjua5JVJ3pfkniTvaq3dVVWvr6qXDIr92yTnJXl3Ve2uqlsWuDsAAFiRRrHGOa21W5PcOmfb62ZdfvEo6uH0MTXdcvueg7nrwOE8Z8uG7Lx8cybWzLdcHgBgZRhJcIbFmJpuecVNH8zufYdyZHIq69dNZPvWjbn5hhcKzwDAiuUnt1l2t+85mN37DuXxyam0JI9PTmX3vkO5fY/l7wDAyiU4s+zuOnA4Ryanjtt2ZHIqdx84PKYWAQCcnODMsnvOlg1Zv27iuG3r103kyi0bxtQiAICTE5xZdjsv35ztWzempiaTNp1zB2ucd16+edxNAwBYkODMsptYU7n5hhdm073/Ixv3/3n+/XXPc2AgALDiOasGYzGxpnLuofty7qH7ctUVF467OQAAJyU4s2TOxQwAnEkEZ5bEuZgBgDONNc4siXMxAwBnGsGZJXEuZgDgTCM4syTOxQwAnGkEZ5bEuZgBgDON4MySOBczAHCmcVYNlsy5mAGAM4kZZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA7OqsGqMzXdcvueg7nrwOE8Z8uG7Lx8s9PgAQCnnODMqjI13fKKmz6Y3fsO5cjkVNYPfnjFOaQBgFPNUg1Wldv3HMzufYfy+ORUWpLHJ6eye9+h3L7n4LibBgCc5gRnVpW7DhzOkcmp47YdmZzK3QcOj6lFAMCZQnBmVXnOlg1Zv27iuG3r103kyi0bxtQiAOBMYY3ziP3FRx8ZdxOW1eEjn0uytMe9lH3PWTuRyy54Su762CeSibU5+6y1ueyCp+SctRNnXN8DwOnmRc86f9xNOCEzzqwqa9ZUXnPNFTnv7t/N+r/7s/zY12/La665ImscGAgAnGKCM6vOmjWVdY/szfoH/jzPv+RpQvMpNj3dcucDj+a379yfOx94NNPTbdxNAoCxsFQDWND0dMsb3ntP9h58LJNHp7Nu7Zo8e/N5ZvkBOCOZcQYWtHvfoew9+FieODqdluSJo9PZe/Cx7N53aNxNA4BlJzgDC7r/kc9k8uj0cdsmj07n/kc+M6YWAcD4CM7Agi49/ylZt/b4l4l1a9fk0vOfMqYWAcD4CM7AgrZv3Zhnbz4vOTqZtOmcPVjjvH3rxnE3DQCWneAMLMjp/wDgCwRn4ISc/g8AZjgd3Wlgerpl975Duf+Rz+TS85+S7Vs3CjcAACMmOK9yzrMLQI8zcZLlTHzMnFqC8yo3+zy7yfHn2X3+JU8bc+sAWAmGnWRZjQF03BNLq7HPODnBeZU70Xl2Bef5eTEDzjTDTLKMO4Au1TgnllZrn3FygvMqd+w8u0/MCs/Os7uwM3HWZbXS1zA6w0yyrNZvNsc5sbRa+4yTE5xXuWPn2b3rY59IJtbm7LPWOs/uCZyJsy6rkb6G0RpmkmW1frM5zoml1dpnnJzT0a1yzrO7OMP8hPTs0N1yfOhmtPQ1jNYwP2a0Wn9BdJw/4LRa+4yTE5xPA86z22+YF7NhQjeLo69htIaZZFmtvyA6zoml1dpnnJzgzBnlTJx1WY3O1L6enm6584FH89t37s+dDzya6ek27iZxGlnqJMu4v9kc5nkxromlcfcZp441zpxRjr2Y/eCPvypT512YV/7Qjd0HnVlPvnxWc18v9aBG67pZyY4F0DyyN8+/5KeXrd7V/LwYV59xagnOnHGW+mI2TOhmcVZrXw/zJu8ofHiycZ9Szpl9mEtwhkUwg7B8VmNfD/Mm7yh8eLJxPS9W80w3p5Y1zgAjMsxBjWfqum4WZ7Wug19qu8f1vHBmHxYykhnnqro6yb9LMpHkLa21X55z+9lJfiPJVyR5JMl3tNbuH0XdACvFMOeNHXZdt6+VF2ec/XWmrYMfpt3jOt7BN0AsZOjgXFUTSd6U5BuS7E9yR1Xd0lq7e1axG5I82lp7dlW9PMm/TvIdw9bN8LzZwugM8yY/zLru1RqoxmWc/XUmroMfpt3jOt7Br/KykFEs1XhBkr2ttftaa5NJ3pnk2jllrk3y9sHl9yS5qqq8mo/ZsRfwN77/3rznQ/vzxvffmze8955V89UfrDTDnoJqqafO8rXy4oyzv4ape7We33zYdo/jlHLjPg/zal2Scyao1ob7z6iqlyW5urX2/YPrr0jywtbaK2eV+etBmf2D6x8dlPnEQvf79EuuaN/wmrcO1bal2P2R3UmS7V++fUn7H/7s50bZnG733v3XSZJtVz63e59Pf/ZoPn7oSGYPgarkmRvX56nn9H0ZsZR6x73vuOpureWxJ6by2c9N5ZyzJnLe2RNZLZ8fh+3v1VbvsJZ7fD386Sfyiccmn7R903nrcsFTz150GxZrtY3tcfbXMHWv1tfs1dru1lr+du99ycS6bNnyjCWN66XU3VrLxz55JEc+N5XWZvpq/VkTufjp61f082pYx15H0pJzz57IxvVnLevjfdcPfdWHWms7TlZuFGuc53tUc9N4T5lU1Y1JbkyS857xrOFbtgRLDczJTOiemmpjCYJL2e+zgyflbK0lT3xuqvvFbJhAM659x1H3sRfCx5/4XJJKralFvxCO88PCMP21Gusddv/lHl/nnDWRqjwpmJx91kTX/sMEhHGP7aXsO2x/javu886eyPqzJp4UqM47u7/d43jdXa3trqpcvm24LLKUuh97YurzfZXMjJUjn5vKY0/0vzevtkmtua8ja9ZUzjt7bf7BFz91xX1YGMWM84uS/Hxr7X8bXP+ZJGmt/atZZd43KPMXVbU2yd8n2dROUPmOHTvarl27hmrbctu5c2cOH/lc3vRbtyx63x/5zpckyZL2Xao7H3g0b3z/vcet4Tp77Zr82NdvW9Hr5VajUfT1sGNkHGNsNdc7rnYvxTDrZo/tO3dddu8Sk3GP7aXsO6o1zuOoe7Uel7Ja2z0Ov33n/rznQ/uPm12sJC/7iovyrc+/qOs+lvs5Ney+872OnLtuIv/+uuflqisuXHQ7lqKqlm3G+Y4k26rqsiQfT/LyJN85p8wtSa5P8hdJXpbk/ScKzSyPY2u45r6Ar4ZfZ1ttHKHNqXRsbfVSgsmxNbdZuy7J4g84W41je5j+OmZ6umXy/Gdn6rwLc+cDj3bvP2zda9ZUnn/J01Zs3y5ktbZ7HIY9MHGpY3Oc5nsdOTI5lbsPHF624Nxr6ODcWjtaVa9M8r7MnI7ura21u6rq9Ul2tdZuSXJTkpuram+ST2YmXDNmo3jzoI8jtDnVlhpMhg2+q3VsDxPkjs0aP3blS5OJtXnj++9d1KyxEMmJDDOpNezYHMYwgX2+15H16yZy5ZYNp6q5SzaS8zi31m5Ncuucba+bdfmzSb5tFHUxWl7Al4fZfVaqYYPvmTi2h52lhxMZ5zdISzVsYJ/7OrJ+3US2b92YnZdvPmVtXio/uQ3LwOw+K9WwwXecY3tcX0mvxuUprC7j+gZpqc+pYQP77NeR6dZy5ZYN2Xn55kyswPdIwRmWidl9VqJRBN9xjO1xfiW9WpencPobZmwO85waxYfJY68jL3rW+V3lx2UUP4ACwCp27A3rW59/0bL9wMSwjpvhqjXL+iMmx2bpz167JpUs+49jwEKGGZvDPKeOBfbZTtcPk2acAVh1xrlcwtIrVqphxuYwz6kz6VgHwRmAVWfcyyUsvWKlWurYHOY5dSZ9mBScAVh1zqQZLlgOozhQ+Ez4MCk4A7DqnEkzXLAcPKf6CM4ArEpnygwXLBfPqZNzVg0AAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDJx2pqdbJs9/do5c8tW584FHMz3dxt0kAE4DgjNwyowjwE5Pt7zhvffksStfmiOXfW3e+P5784b33iM8AzA0wRnOAGdSgN2971D2HnwsWbsuqTV54uh09h58LLv3HTql9QJw+hOcYRUYJvieaQH2/kc+k8mj08dtmzw6nfsf+cwprZelsawGWE0EZ1jhhg2+Z1qAvfT8p2Td2uNf2tatXZNLz3/KKa2XxbOsBlhtBGdY4YYNvmdagN2+dWOevfm8nL12TSrJ2WvX5Nmbz8v2rRtPab0snmU1wGqzdtwNON1sWH9WXvSs85e0X5Il7cvp7Y77Pzlv8J1urWu8PD55NH/wVw/m8cmpz29bv24iVz/3i0/peHvBZU/Pn3/0E9m971COTE5l/bqJbN+6MT+081mZWFOnrN4kueVZX5Pb9xzM3QcO58otG7Lz8s2LrtNz8tQbdmwDLDfBGVa452zZkPXrJp4UfK/csqFr/52Xb872rRufFGB3Xr75VDU5STKxpnLzDS8cOsAute6rrrgwV11x4Smvi6UbdmwDLDfBGVa4YYOvAMtKNa4PdQBLJTivAFPTLY9v/JJMPuXC3HbPQ8sWalgdRhF8BVhWonF+qANYCsF5zKamW15x0wfz8LZvSVuzNj/6jg9n+9aNufmGF3rz4PME3+Xlw+zyMbaB1cRZNcbs9j0Hs3vfobSJmaPKH5+cyu59h3L7noPjbhqckWZ/mD100VflR9/x4bzipg9myinSAM54gvOY3XXgcI7MOjAmSY5MTuXuA4fH1CI4s/kwC8BCBOcxO3ZU+WyOKofx8WEWgIUIzmN27Kjyc9dNpJKc66hyGCsfZgFYiIMDx8xR5bCyOEUaAAsRnFcAR5XDyuHDLAALEZwB5vBhFoD5WOMMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADr6r/5HAAAJ2UlEQVQIzgAA0EFwBgCADoIzAAB0GCo4V9XTq+qPqurewd+nzVNme1X9RVXdVVV/WVXfMUydAAAwDsPOOL86yW2ttW1Jbhtcn+vxJN/TWntOkquT/HpVbRyyXgAAWFbDBudrk7x9cPntSV46t0Br7W9ba/cOLh9IcjDJpiHrBQCAZTVscL6wtfZgkgz+bj5R4ap6QZJ1ST66wO03VtWuqtr18MMPD9k0AAAYnbUnK1BV/zPJF89z088upqKqekaSm5Nc31qbnq9Ma+3NSd6cJDt27GiLuX8AADiVThqcW2svXui2qnqoqp7RWntwEIwPLlBuQ5I/SPLa1toHltxaAAAYk2GXatyS5PrB5euT/N7cAlW1LsnvJPmN1tq7h6wPAADGYtjg/MtJvqGq7k3yDYPrqaodVfWWQZlvT/KPknxvVe0e/Ns+ZL0AALCsTrpU40Raa48kuWqe7buSfP/g8m8m+c1h6gEAgHHzy4EAANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAh6GCc1U9var+qKruHfx92gnKbqiqj1fVfximzpVqarrl8Y1fkkPPfFFuu+ehTE23cTcJAIARGnbG+dVJbmutbUty2+D6Qn4xyZ8MWd+KNDXd8oqbPpiHt31LDl30VfnRd3w4r7jpg8IzAMBpZNjgfG2Stw8uvz3JS+crVFVfkeTCJH84ZH0r0u17Dmb3vkNpE+uSWpPHJ6eye9+h3L7n4LibBgDAiAwbnC9srT2YJIO/m+cWqKo1SX4lyU8OWdeKddeBwzkyOXXctiOTU7n7wOExtQgAgFFbe7ICVfU/k3zxPDf9bGcdP5zk1tbavqo6WV03JrkxSS6++OLOux+/52zZkPXrJvL4rPC8ft1ErtyyYYytAgBglE4anFtrL17otqp6qKqe0Vp7sKqekWS+tQkvSvK1VfXDSc5Lsq6qHmutPWk9dGvtzUnenCQ7duxYNQuEd16+Odu3bszufYdyZHIq69dNZPvWjdl5+ZMm4AEAWKVOGpxP4pYk1yf55cHf35tboLX2XccuV9X3JtkxX2hezSbWVG6+4YW5fc/B3H3gcK7csiE7L9+ciTUnnmEHAGD1GDY4/3KSd1XVDUk+luTbkqSqdiT5odba9w95/6vGxJrKVVdcmKuuuHDcTQEA4BSo1lbmiogdO3a0Xbt2jbsZAACc5qrqQ621HScr55cDAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0KFaa+Nuw7yq6uEkD4yp+guSfGJMda9G+mtx9Nfi6K/F0V+Lo78WR38tnj5bnHH11yWttU0nK7Rig/M4VdWu1tqOcbdjtdBfi6O/Fkd/LY7+Whz9tTj6a/H02eKs9P6yVAMAADoIzgAA0EFwnt+bx92AVUZ/LY7+Whz9tTj6a3H01+Lor8XTZ4uzovvLGmcAAOhgxhkAADoIzrNU1dVVtaeq9lbVq8fdntWgqu6vqr+qqt1VtWvc7VlpquqtVXWwqv561ranV9UfVdW9g79PG2cbV5IF+uvnq+rjgzG2u6q+cZxtXEmqamtV/XFV3VNVd1XVjw+2G2PzOEF/GWPzqKpzqur/q6qPDPrrFwbbL6uqDw7G13+rqnXjbutKcIL+eltV/d2s8bV93G1dSapqoqo+XFW/P7i+oseX4DxQVRNJ3pTkmiRXJrmuqq4cb6tWja9rrW1fyaePGaO3Jbl6zrZXJ7mttbYtyW2D68x4W57cX0nya4Mxtr21dusyt2klO5rkVa21K5J8ZZIfGbxuGWPzW6i/EmNsPk8k+frW2pcn2Z7k6qr6yiT/OjP9tS3Jo0luGGMbV5KF+itJfnLW+No9viauSD+e5J5Z11f0+BKcv+AFSfa21u5rrU0meWeSa8fcJla51tqfJvnknM3XJnn74PLbk7x0WRu1gi3QXyygtfZga+3OweVPZ+bN55kxxuZ1gv5iHm3GY4OrZw3+tSRfn+Q9g+3G18AJ+osFVNVFSb4pyVsG1ysrfHwJzl/wzCT7Zl3fHy+oPVqSP6yqD1XVjeNuzCpxYWvtwWTmjTzJ5jG3ZzV4ZVX95WAph2UH86iqS5M8L8kHY4yd1Jz+SoyxeQ2+Rt+d5GCSP0ry0SSHWmtHB0W8V84yt79aa8fG1y8NxtevVdXZY2ziSvPrSX4qyfTg+vlZ4eNLcP6CmmebT4on99WttednZonLj1TVPxp3gzjt/Mckz8rMV58PJvmV8TZn5amq85L89yT/vLV2eNztWenm6S9jbAGttanW2vYkF2Xmm9kr5iu2vK1aueb2V1U9N8nPJPkHSf5hkqcn+ekxNnHFqKpvTnKwtfah2ZvnKbqixpfg/AX7k2yddf2iJAfG1JZVo7V2YPD3YJLfycwLKyf2UFU9I0kGfw+OuT0rWmvtocGb0XSS/xxj7DhVdVZmQuB/ba399mCzMbaA+frLGDu51tqhJLdnZm34xqpaO7jJe+U8ZvXX1YMlQq219kSS/xLj65ivTvKSqro/M8tjvz4zM9ArenwJzl9wR5Jtg6M51yV5eZJbxtymFa2qnlJVTz12Ock/SfLXJ96LzIyr6weXr0/ye2Nsy4p3LAAO/NMYY583WA94U5J7Wmu/OusmY2weC/WXMTa/qtpUVRsHl9cneXFm1oX/cZKXDYoZXwML9NffzPoQW5lZr2t8JWmt/Uxr7aLW2qWZyVzvb619V1b4+PIDKLMMTkH060kmkry1tfZLY27SilZVX5KZWeYkWZvkt/TZ8arqHUl2JrkgyUNJfi7J7yZ5V5KLk3wsybe11hwQlwX7a2dmvkJvSe5P8oPH1u+e6arqa5L8WZK/yhfWCL4mM+t2jbE5TtBf18UYe5Kq+rLMHJw1kZmJtne11l4/eO1/Z2aWHXw4yXcPZlPPaCfor/cn2ZSZZQi7k/zQrIMISVJVO5P8H621b17p40twBgCADpZqAABAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOjw/wOCsl1szxo+SgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 864x576 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "#使用一阶差分，12步差分处理时间序列\n",
    "diff_1 = data.diff(1)\n",
    "diff1 = diff_1.dropna()\n",
    "diff1_144_1 = diff_1-diff_1.shift(144)\n",
    "diff1_144 = diff1_144_1.dropna()\n",
    "#print(diff1_144_1)\n",
    "#判断序列是否平稳，计算ACF，PACF\n",
    "fig1 = plt.figure(figsize=(12,8))\n",
    "ax1=fig1.add_subplot(111)\n",
    "sm.graphics.tsa.plot_acf(diff1_144,lags=40,ax=ax1)\n",
    "fig2 = plt.figure(figsize=(12,8))\n",
    "ax2=fig2.add_subplot(111)\n",
    "sm.graphics.tsa.plot_pacf(diff1_144,lags=40, ax=ax2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#模型定阶，根据aic,bic,hqic,三者都是越小越好\n",
    "# arma_mod01 = sm.tsa.ARMA(diff1_144,(0,1)).fit()\n",
    "# print(arma_mod01.aic,arma_mod01.bic,arma_mod01.hqic)\n",
    "# arma_mod10 = sm.tsa.ARMA(diff1_144,(1,0)).fit()\n",
    "# print(arma_mod10.aic,arma_mod10.bic,arma_mod10.hqic)\n",
    "# arma_mod60 = sm.tsa.ARMA(diff1_144,(6,0)).fit()\n",
    "# print(arma_mod60.aic,arma_mod60.bic,arma_mod60.hqic)\n",
    "arma_mod61 = sm.tsa.ARMA(diff1_144,(6,1)).fit()\n",
    "print(arma_mod61.aic,arma_mod61.bic,arma_mod61.hqic)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.5.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
